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Abstract. The survey provides classification results for integrable one-field 
evolution equations of orders 2, 3 and 5 with the constant separant. The classifi- 
cation is based on necessary integrability conditions following from the existence 
of the formal recursion operator for integrable equations. Recurrent formulas 
' for the whole infinite sequence of necessary conditions are presented for the first 

Q ■ time. The most of the classification statements can be found in papers by S.I. 

CN . Svinilupov and V.V. Sokolov but the proofs have never been published before. 

^ i The result concerning the fifth order equations is stronger than obtained before. 

Keywords: evolution differential equation, integrability, generalized symmetry, 
conservation law, classification. 
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Introduction 

This survey is devoted to the classification of integrable evolution equations 

(0.1) Ut = Un + F{X,U, U^, Uxx, ■ ■ ■ ,Un-l>, ^ 

The equations with such dependence of the highest x-derivative are often referred to as equations 
with a constant separant. 

Let us specify what we mean by the integrability in the present paper. Unfortunately, at 
present there exists no unified rigorous definition for the integrability of differential equations 
(for various approaches see, for instance, [IH3]). However, for some types of differential equa- 
tions there are efficient criteria of the integrability, which can be not only checked for these 
equations, but also allow one to find all the equations from this class satisfying this criterion. 

For evolution equations (10. ip with one temporary variable and one spatial variable the most 
effective integrability criterion is the existence of generalized local symmetries. In the works [IfT] 
a way of "excluding a symmetry" from this relation and obtaining necessary conditions for the 
existence of symmetries only in terms of the right hand side of the equation was suggested. 
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These conditions which we call integrability conditions are written as so-called canonical con- 
servation laws. Their main advantages are the independence of the conditions on the order of 
symmetry and their invariancy w.r.t. all point transformations not leading out of the class of 
equations flO.ip . 

It was shown in the papers [1H7] that necessary integrability conditions are implied by the 
existence of an infinite series of generalized symmetries or conservation laws for equation flO.ip . 
In more details the technique of obtaining the conditions is discussed in the reviews [SHI]- 
Here we do not deal with it. We note that there exists an alternative way [IOl[TT] for calcu- 
lating canonical conservation laws by the logarithmic derivative of the formal eigenfunction of 
the linearization operator for equation (10. ip (see Appendix 3). The equivalency of these two 
approaches for the scalar equations follows from Theorem 2.9 in survey [T^ . 

Let us describe the results of the work. In Chapter 1 by the simplest examples we show 
how canonical conservation laws look like and how one can classify integrable equations by 
employing them. In particular, in this chapter the problem of the classification for equations 
(10. ip with n = 2 is solved. General second order integrable evolution equations were classified 
in [13] . In [13] the results of the latter work were generalized for the case of weakly nonlocal 
symmetries. 

In Chapter 2 we provide the solution of the classification problem for the integrable equations 
of the form 

(0.2) ut = us + F{x,u, ui, U2). 

The famous Korteweg-de Vries equation 

(0.3) Ut = U2, + uui 

belongs to this class. The case when the function F is independent on U2 and x (see Section 
1.2) was considered in [HHS]. The results of Chapter 2 were announced in [SIE], but the proof 
is published now for the first time. We also present for the first time a recurrent formula 
describing all infinite series of canonical densities. In the works OE] only 4 first densities were 
written down explicitly which were indeed used then in the classification. Third order evolution 
integrable equations more general than (10. 2p were studied in [9t[T6|fT7]. 

In Chapter 3 we consider a computationally complicated problem on the classification of 
integrable equations of the form 

(0.4) = ^5 + F(^^,^^l,^^2,^^3,^^4)• 

In the note [18j a solution to this problem was announced under an additional assumptions 
that even canonical densities are trivial (see Remark 2). However, not only the proof but 
also any complete list of the found equations is absent in [IB]. For the first time the list of 
equations (10.41) possessing generalized conservation laws was published in [9]. In the present 
work the condition of the triviality of even canonical densities is not employed and we solve 
thus a technically more complicated problem on the classification of equations (10.41) possessing 
generalized symmetries. The answer coincides in essence with the list in [D]. As in the case of 
the third order equations, a general formula for the whole infinite series of canonical densities 
is published for the first time in the present paper. 

The results of the works [5|[^[^[TB] were obtained by hard calculations made "by hand" . This 
is why it was a non-zero probability of errors which could lead to losing integrable equations. 
Once computer systems like Maple, Mathematica, etc. appeared, an opportunity to automate 
partially the calculations rose. The results of the present paper were obtained by the program 
package Jet written by the first author. It was found no essential errors in the lists of the 
integrable equations but we found and corrected several misprints in [9]. 
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At first glance, the problem of the classification of integrable equations (10. ip with arbitrary 
n seems to be far from the complete solution. This is not quite so. Each integrable equation 
together with all its symmetries form a so-called hierarchy of integrable equations. For the 
equations integrable by the inverse scattering problem method [12] all the equations of the 
hierarchy possess the same L-operator. This fact lies in the basis of the commutativity of the 
flows in hierarchies (each equation of the hierarchy is a symmetry for all others). A general 
statement on "almost" commutativity of the symmetries for equation (10. ip is contained in [2U] . 

Assuming that the right hand side of equation (10. ip is polynomial and homogenous, it was 
proven in the works [211|22] that the hierarchy of any such equation contains an equation of 
second, third, or fifth order. This statement looks very credible also without any additional 
restrictions for the right hand side of the equation. The proof in the general case is absent and 
this statement has a status of the conjecture well-known for experts. No counterexamples to 
this conjecture are known. 

Up to this conjecture, in the survey we describe all the hierarchies of the integrable equations 
of the form (10. ip . In other words, any integrable equation of order 4 or > 5 is equivalent to a 
generalized symmetry of one of the equations given in this survey. We note that the calculation 
of symmetries for given equation is a linear problem, and there are several effective computer 
programs for solving it. Moreover, the generalized symmetries can be found by the use of 
quasilocal recursion operators (see |23] and the references therein). 

Various results on the classification of integrable systems of evolution equations can be found 
in [5| I21] - HD] . Further references are contained, for instance, in the survey [H]. 

A separate difficult problem is the classification of integrable hyperbolic equations and sys- 
tems [12H50]. 

Acknowledgments. The authors are grateful to D.I. Borisov, A.V. Mikhailov, S.I. Svinolupov, 
and A.B. Shabat for numerous useful discussions. V.S. is grateful to Max Planck Institute 
(Bonn) for the hospitality. The research is partially supported by RFBR grant 11-01-00341- 
a, the grant for supporting leading scientific school 6501.2010.2, and the grant of Ministry of 
Education and Science of Russia (project 1.2.11). 



1. Simplest classification problems 

All necessary integrability conditions we shall use below are given in the form of local con- 
servation laws. We remind [HT] that a pair of functions p and 9 depending on a finite number 
of the variables x,u,Ui, . . . such that 

is called a local conservation law for equation (10. ip . Here 

d d d d d 

(1-2) 1- = ^+"!^ ^"2^ hM3^ \ , Uq = U, 

ax ox ouq oui 0U2 

d d d d d 

-77 = it: + Ko— + Ki— + R2- 



dt dt duQ dui du2 
where 



d^ f \ 
-f^j = ^ V " ■ ■ ■ ' ^"-1) j • 
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the operators and are often referred to as the total derivative w.r.t. x and the total 
dx dt 

derivative w.r.t. t in virtue of equation (10. ip . The function p is called a density, and 6 a flux 
of the conservation law. 

Relation (11. ip is called a conservation law due to the following reason. Consider, for instance, 
the Korteweg-de Vries equation ut = + uui. It is known that it possesses an infinite number 
of conservation laws. In particular, since the equation can be rewritten as 

Ut = {U2 + 

the function u is the density of the conservation law. Suppose the solution u{x,t) decays as 
|x| — )■ oo. Then we have 

d /•+~ , 

— / udx = 0, 

dt i_oo 

i.e., the area under the graph of the solution is independent of t. In the same way, the integrals 
of others densities of the conservation laws are conserved. ^ 

It is clear that if p is a density of a conservation law, then pi = p + — {h) is also a density 

dx 

for any function h. We call two such densities equivalent and write p ~ pi. A conservation law 
is called trivial if p ~ 0. 

The order of the higher derivative, on which the function /(x, m, mi, . . . , Uk) depends, is called 
differential order of this function. The differential order is usually indicated as ord f = k. The 
minimum of the differential orders of equivalent densities is called the order of the conservation 
law. 

The deduction of necessary integrability conditions as an infinite series of so-called canonical 
conservation laws was discussed in details in [5HTT]: for an alternative version see Appendix 3. 
In this paper we often give appropriate formulas without proofs. But on the other hand we 
dwell on how to retrieve the complete list of integrable equations of the form (10. 2p using these 
necessary conditions, and we describe point transformations necessary to reduce an arbitrary 
integrable equation to one of the canonical forms. 

1.1. Integrable Burgers type equations. Consider second order evolution equations 

(1.3) Ut = U2 + f{x, U, Ui). 

The canonical densities for this equation are defined by the recurrent formula 

(1.4) 2p„+i = 6^„, + 2_^p„_iPi - — p„ + — + — d„o - ^Pn, 

i=0 ^ ^ 

Here p_i = 0, 6ij is the Kronecker delta. One of the ways of obtaining similar formulas is 
described in Appendix 3. The fluxes associated with these densities are calculated consequently 
in the process of classification. At that, the obstacles to their existence pose the restrictions for 
the right hand side of equation (10. 2p that finally allow us to find all integrable equations (II. 3p . 
Letting n = —1,0 in (II. 4p . we find two first canonical conservation laws, 

d df d 

(1.5) 

ui uui ax 

d ( df 

2 \dui J I dx 



where ai = 2 6*0 and a2 




1.6) i-M^-r^.-Um-i"^ 
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The former of these formulas means that for each integrable equation (11. 3p the partial de- 
rivative w.r.t. Ui of its right hand side is a density of a conservation law. For instance, for 
Burgers equation Ut = U2 + uui this formula yields a density p = u. In this case the function 
0"! is calculated easily, 

1 2 

(yi = u2 + -u . 

A general algorithm of calculating the flux for a given density is given below (see Remark 4). 

Let us demonstrate the main modes for working with the conditions like (II. 5p . (II. 6p . In order 
to determine the character of the dependence of the right hand side on ui, the most simplest 
way is to exclude the unknown function o"i in (II. 5p . For this we apply the Euler operator 

5 d d d (f d 
o 1 o 



5u du dx dui dx^ du2 
to both sides of (II. 5p . It is well known [ST] that 



6 d 
6u dx 



0, 



and therefore, 

(17) = -- (^] = -2uJ^ - 4^3-^ + 0(2) 

5udt\duiJ du\ dxdu\ 

where the symbol 0(2) indicates terms whose order w.r.t. the derivatives are at most two. The 
latter identity must hold true for each solution (11.31) . Since there exists no ordinary differential 
equation in x satisfied by all the solutions to equation (II. Sp . relation (ll.7p must hold identically 
w.r.t. the variables u,ui, . . . ,U4. Equating the coefficient at M4 to zero, we find that the 
equation reads as 

(1.8) Ut = U2 + A{x, u)ul + B{x, u)ui + 0(x, u). 

Thus, each integrable equation (11.31) is quadratic in ui. It can be checked that for equation 
(ll.Sp condition (11.71) is equivalent to two equations 

where ip = Bu — '^A^. 

Taking into consideration that the integrability of any differential equation is preserved under 
point transformations, we simplify equation (II. 8p by a point transformation u = ip{x,v) before 
we proceed to study the integrability conditions. Simple calculations lead us to the evolution 
equation for v 



It is obvious that the equation 

has a solution depending on v forany function A. This is why by a point transformation one 
can vanish the function A in equation (II. 8p . This transformation is the first step in reducing 
any integrable equation to one of the canonical forms. 
Condition (II. 5p for the equation 

(1.9) Ut = U2 + B{x,u)ui + C{x,u) 
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becomes 

(1.10) Bu{u2 + B{x,u)ui + C{x,u)) = -^CTi- 

Since to use condition f ll.6p we need to know completely or partially the function ui, instead 
of applying the variational derivative to both sides of fll.lOp we employ an alternative approach 
which is a separation of a total x-derivative in the left hand side of fll.lOp . This approach is 
completely algortihmical and can be programmed in any language of symbolic computations 
(see Remark on page ITT]) . 

We have 

5„M2 + B^Bui + BuC = ^ (^B^ui + ^B^^ - Mi(5„„Mi + 5^ J - BB^ + B^C = 
— ( BuUi + -B^ — Bx \ — Buu'ul + Bxx — BBx + -B«C. 



dx\ 2 
Substituting the last expression in fll.lOp . we obtain 

~B^uu\ + 5... - BB^ + B^C = ^ (^1 - ^."1 - + B.)j = ^. 

Since the left hand side depends only on x, u, ui, then the function if) may depend only on x 

and u. Substituting — ^ = tpx + i'uUi and equating the coefficients at u\ and ui, we get 

dx 

Buu = 0, V'n = 0, Bxx — BBx + BuC = ipx- 
Letting B = a{x)u + (3{x), we find that each integrable equation (11. 9p reads as 

(1.11) Ut = U2 + {a{x)u + /3(x))mi + C{x, u), 
where 

(1.12) aC{x, u) - aa'u^ + {a" - a/3' - a' P)u = tp' + (3(3' - (3". 
At that, 

(1.13) ai = ip + aui + -{au + (3)'^ — a'u — (3'. 

If a 7^ 0, then by fll.l2p one can determine the function C. In this case equation (11. lip can be 
simplified by a point transformation u u fi{x)+f2{x). By taking fi = 1/ a = ^a' /a^ — (3 /a 
we get a = 1, /9 = 0. At that, equation (II. lip casts into the form 

(1.14) Ut = u^x + uu,j. + ^()'{x). 

For this equation conditions (II. 5p . (II. 6p . as well as all other necessary integrability conditions, 
hold true. The Burgers equation (I1.14p is reduced to the linear equation 

by the Cole-Hopf substitution u = 2vx/v + f{x), where ip and ip are related by the identity 

^" + = - f. 

In the case a = the left hand side of equation (ll.lOp vanishes and this is why ai is constant. 
Then condition (II. 6p implies 

= 0, 





INTEGRABLE EVOLUTION EQUATIONS WITH CONSTANT SEPARANT 7 

that is equivalent to the system of equations 

Cuuu = 0, CCuu + Cxuu ~ {/^Cu)x + ^ {x) = 0. 

It yields C = p{x)u + q{x), and we arrive at the hnear equation 

(1.15) Ut = + I3{x)ux + p{x)u + q{x). 

For this equation all necessary integrability conditions hold true. 

Remark 1. Among obtained second order integrable equations (11.141) and (ll.lSp there is 
no the potential Burgers equation Ut = u^x + u"^- The reason is that this equation is linearized 
by the point transformataion m = In^;. This transformation is a particular case of the point 
transformation that has to be applied to equation (II. 8p for eliminating the function A. 

1.2. Integrable KdV type equations. The list of integrable equations obtained in the pre- 
vious section is quite poor. Let us consider a more substantial classification problem. Let us 
find all integrable evolution equations of the form 

(1.16) Ut = us + f{ui,u). 

It turns out (see Section 2.1) that for each such integrable equation 

where ai is a function depending on u, Ux, . . . , u^. 

Example 1. For the mKdV equation ut = U3 + u^ui conservation law (I1.17P reads as 

{u\ = {2uu2 -u\ + \l^)x. □ 
Applying the Euler operator to both sides of (I1.17p . we obtain 

The last identity must hold true for each solution of (I1.16P and this is why it must be identity 
in the variables . . . ,^4. Equating the coefficient at M4 to zero and employing that / is 
independent of u^, we get 

f{ui,u) = fiul + A{u)uj + B{u)ui + C{u) 

with some constant /i. It is easy to check that for each such function / condition (11.181) is 
equivalent to the system of ODEs 

/iA' = 0, 5"' + 8/i5' = 0, (5'C)' = 0, A5' + 6/iC" = 0. 

The next necessary integrability condition reads as 

d / df\ d . . 

that yields 

^ ^ Su dt ^dii^ 

The latter condition leads to additional equations 

A' = 0, AC" = 0, (C" + 2/iC")' = 0, {CC")' = 0. 
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In the case fi ^ the obtained equations are sufficient to determine completely the functions 
A, B, and C. As a result, up to a scaling u — )■ const u, we arrive at the equations 

(1.20) ut = M^^^ - ^ul + (cie^" + c2e-2" + Cg) 
and 

(1.21) Ut = U^xx + Clul + C2ul + CsM^ + C4, 

where Cj are arbitrary constants. 

If /i = 0, then solving the above system of ODEs for the functions A, B, C, we obtain that 
the equation reads as 

Ut = Uxxx + CquI + {ciU^ + C2U + C3)Ux + C^U + C5, 

where 

CqCi = 0, C0C2 = 0, C4C1 = 0, C4C2 = 0, C1C5 = 0. 
By the third integrability condition (see Section 2) we find additional relations, 

C0C4 = 0, C2C5 = 0. 

In the case cq 7^ we arrive at a particular case of equation fll.211) . If cq = 0, two cases are 
possible; a) ci 7^ or C2 7^ 0, C4 = C5 = and b) ci = C2 = 0. They lead us to two equations 

(1.22) Ut = U^xx + (ciM^ + C2U + 03)^0:, 

(1.23) Ut = Mxxx + CsMx + C4U + C5. 

The experts in nonlinear equation regard each linear equation as exactly integrable. Equa- 
tions fll.20p . (11.211) . and (11.221) have been found by necessary integrability conditions. This is 
why it should be discussed independently in which exactly sense they are integrable. It is well- 
known that to all of these equations the method of the inverse scattering problem is applicable. 
Moreover, all of them are related with the Korteweg-de Vries equation Ut = U3 + uui by Miura 
type differential substitutions [52] . 

Remark 2. Conditions (11.181) . (11.191) hold true for the equations (I1.16P possessing generalized 
symmetries. If the equation possesses generalized conservation laws (at that, the existence of the 
symmetries is not assumed), condition (11.181) still holds, and condition (I1.19P can be strengthen, 

5u \du 

It is implied by the general statement |5] in accordance to which for the equations with gener- 
alized conservation laws the canonical densities with even indices are trivial. 

1.3. On admissible point transformations. In the process of classification of integrable 
equations, as a rule, we use point transformations reducing integrable equation to one or another 
canonical form. For instance, in Section 1.1 we employed point transformations while reducing 
equation (11.80 to the form (II. 9p . and also for normalizing the functions a{x) and /3(x) in 
equation (11.111) . 

Let us describe point transformations we use in the classification of equations (10. ip . 
Each equation of the form (10. ip admits the transformations 

(1.24) u = ip{u,x). 

Hereinafter, once transformation rules for some of the variables t, x, or u are not indicated in 
the formulas, this means that the corresponding variables are not changed. The scalings 

(1.25) x = ax, i= aJ't 
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are also admitted. Under such transformations 

F{x, u, Ui, U2, . . .) — )■ a~"'F{a~^x, u, aui, a^U2, . . . ). 

For some subclasses of equations ( 10. ip additional transformations depending on t are ad- 
mitted. In particular, if Am, Ami, • • • , Am„_i) = XF{x, u,ui, . . . , Un-i), then for arbitrary 
constants a and b the transformation 

(1.26) u = uexp{at + bx) 

is applicable. Under this transformation Un {dx — b)"'u, F ^ F + au. 

If, as in Section 1.2, it is assumed that the right hand side F of equation flO.ip is independent 
on the variable x, then the class of admissible transformations changes. Among fll.24p . only 
the transformations 

(1.27) u = ip{u) 

are admitted. At the same time additional point transformations appear. In particular, the 
Galilean transformation 

(1.28) x = x + ct 

is always admissible; under this transformation F ^ F — cui. li the function F is independent 
of u and x, then the transformation 

(1.29) u = u + CiX + C2t 
is admissible. Under such transformation 

F{ui, U2, Ms, . . . ) F{ui - Ci, U2, Ms, . . . ) + Cg. 

The equations related by aforementioned transformations are called equivalent. It is impor- 
tant to note that our classification is pure algebraic. We are not interesting in such properties 
of the solutions to the studied equations as being real. This is why the functions and con- 
stants being involved in formulas fll.24p -( ll.26p can be both real and complex. For instance, 
the equations ut = — u\ and Mt = Ms + m^ are regarded as equivalent. 

Integrable equations can involve arbitrary constants which can be eliminated by one or an- 
other transformation. Consider as an example equation fll.2ip . where Ci 7^ 0. By (possible com- 
plex) scaling m —t- Am we fix a normalization ci = 1. Then the transformation m — )■ m + ax + fit 
leads us to the equation 

Mt + /3 = Uxxx + {Ux + af + C2{Ux + of + C3(m^ + tt) + C4. 

It is easy to see that taking a = — C2/3 and (3 = C4 + + C20? + csa, we obtain C2 = 0, C4 = 0. 
The constant Cs is eliminated by the Galilean transformation, and we obtain the potential 
modified Korteweg-de Vries equation, 

s 

^xxx ~l~ ^x' 

Similarly, the parameters in equation (ll.22p are inessential. 



2. Third order equations with constant separant 

2.1. Integrability conditions. For the equations of the form (10. 2 p an infinite chain of the 
canonical conservation laws 

(2.1) l^^'^) = ^ = 0' 
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can be defined by the formulas (for the deduction see Appendix 3), 

d 



Pn+2 



1 " 



n+l 



Ps Pk Pn- 



d 

dx 



Pn+l 



(2.2) -i E 

where the first two elements of the sequence pi read as 



1 d 

^^PsPn-s + --^Pn 

s=0 



s=0 



n ^ 0, 



If 

3 



Pi 



9 



1 1 
3 "^ + 3rf^ 



Here 5i^j is the Kronecker delta, F^j. = dF/dui, where i = 0, 1,2. The fluxes On are calculated 
consequently in the process of classification. At that, the obstacles for its existence lead to 
differential equations, which must be satisfied by the right hand of integrable equation (10. 2p . 
It is easy to check that first four conditions in this series are equivalent to the conditions 



(2.3) 
(2.4) 

(2.5) 
(2.6) 



d OF 
dt duo 



d 

dx 



d_ 

dt 

d_ 
dt 

d_ __ 

dt ^ dx ^ 



dui 

9ao + 2 
d 



dF 



dF 



d 



c^i, 



dx 

dF 
duo 



OF 
dui 



+ 27 



dF 
du 



d_ 

dx 



d 



given in |6|. At that, do = — 36^0, oi = 3— (Tq — 96i, ... As it will be shown below, these four 

dx 

conditions are "almost" sufficient to obtain the complete list of integrable equations (10. 2p . 

In order to employ efficiently the canonical series for the classification, it is useful to study 
first a possible structure of the densities of local conservation laws of small orders for the 
considered class of equations. 

Lemma 1. If a density p of a conservation law for equation liO.S^) has the differential order 
ord p = 2, then 

(2.7) p = fml + /2M2 + /s, 

where fi are some functions in x,u, Ui, and 

(2.8) 



dx 3 0U2 



Proof. Eliminating the terms U5 and U4 by the subtraction of total x- derivatives, we find 
that 

d' 



d 
~dt 



Ui + —F 
dx 



dp 



(2.9) 



u. 



d^p 



2 dui 2 



Un 



d^p d^p 



U5 + 
d'^p 2 dF d'^p 



duldui 



duldu 



duidx 3 du2 dui 



+ 
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where the dots indicate a hnear in U3 expression. By the definition of the conservation law, 

the last expression should read —cr. It is clear that the function a can not depend on the 

ax ^ 

derivatives higher than ti2, and the degree in M3 of the function —a is at most one. Hence, 

d x 

equating the coefficients at u\ and u\ to zero, we obtain fl2.7p and f l2.8p . □ 

Remark 3. We observe that in fl2.7p the identity /i = is possible. It concerns also other 
similar lemmas. 

Remark 4. In the proof of Lemma 1 we used the following algorithm of checking whether 
a given function S{x, u,Ui, . . . , Un) is a complete derivative w.r.t. x (i.e., whether it belongs to 

Im— ). At first, S must be linear in the highest derivative m„. If it holds true, then as one 

iXiJu 

can see easily, we can subtract a total derivative from S so that the difference has the order 
less than n. Repeating this order lowering procedure, we either arrive at the situation when 
the function is nonlinear in its highest derivative, or we get zero. 

Let us show how one can employ formulas (12.70 and (12. 8p in the classification of equations 

Lemma 2. Let for equation (10.21) the first integrability condition (12.31) holds. Then F is a 
polynomial in U2 of at most second degree. 

dF 

Proof. According to condition f l2.3p . the function — — should be a density of a conservation 
law. Applying Lemma 1 to it, we write equations (12. 7p and (12. Sp . 

dF 

= fiul + f2U2 + /s, 

dfi dfi dfi 22,. , . N 

OX OUq OUi 3 

Since fi are independent of U2, then equating the coefficients at u^, we obtain /i = 0. Integrating 
dF 

the equation 7- — = /2M2 + /s w.r.t. U2, we arrive at the desired result. □ 
0U2 

2.2. List of integrable equations. Our main aim is to prove the following statement [6]. 

Theorem 1. Up to transformations of the form (ll.24p - (ll.29p each equation ()U.2p satisfying 
integrability conditions (12. ip . (12. 2p for n = 0, 1, 5, belongs to the list 



(2.10) 


Ut 


U'xxx 


+ UUx, 


(2.11) 


Ut 


Uxxx 


+ U^Ux, 


(2.12) 


Ut 


Uxxx 




(2.13) 


Ut 


Uxxx 


- Kil + (cie'" + C2e-2")M„ 


(2.14) 


Ut 


Uxxx 




(2.15) 


Ut 


Uxxx 


- 2u, + - ^^^^ 


(2.16) 


Ut 


Uxxx 




(2.17) 


Ut 


Uxxx 


2Ux ' 
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(2.18) = Wxxx - ^ + + C2U^, ci 7^ or C2 7^ 0, 

(2.19) ut = u^^^ - + a{x)u^, 



(2.20) 



Q 3/2 



+ (e«" - 2a'2) + / + Co + ci M + C2 



1 . , . ^ a'" , 3a"2 , „a" ^a'^ cq + ci a + C2 
where ^ = a(x) - u, / = y + -—r + 3— - 3— , 

(2.21) «,^,.„ + 3A...9.„^ + 3A. + «.aM + iwW. 

(2.22) Ut = u^xx + 3 UU.J..J. + 3 u^. + 3 u^u^ + {u-f{x))x + f3{x), 

(2.23) Ut = u^xx + a(x)u^. + I3{x)u. 

Here {p'Y = ip^ — g2p — gs, ai, a2, Cq, Ci, C2, g2, 93 are arbitrary constants, a, (3, and 7 are 
arbitrary functions. 

Remark 5. Quite often instead of equations (12.151) and fl2.16p one considers point equivalent 
to them equations 

(2.24) Ut = w... - I ^ + 
and 

(2.25) ^* = ^--^M^T^ + ^^""- 

In both cases, Q = Cq + Ciu + C2U^ + c^u^ + c^u^ is an arbitrary polynomial. If Q' 7^ 0, then one 
can make the substitution u = f{y) in equations fl2.24p and (12. 25 p . where {f'Y = Q{f)- Then 
for V we get equations (I2.15P and (I2.16p . respectively. At that, 

4 2, 8 o 4 32 .2.2 

92 = -C2 - 4ciC3 + I6C0C4, 93 = - -C1C2C3 - y C0C2C4 + 4:CoC^ + Ac^d. 

We observe that under the linear fractional transformations 

(2.26) u = — z 

Z3U + ^4 

the polynomial Q changes by the law 

Q{u) = Q I ) (^3" + ^4)^(21^4 - ^22:3)"^- 

The expression 92,93 are invariants of transformations group (I2.26p . Subject to the structure 
of the multiple roots, by a transformation (12.260 and scalings of x and t the polynomial Q can 
be reduced to one of following canonical forms, Q{x) = x{x — l)(.x — k), Q{x) = x{x — 1), 
Q{x) = x^, Q{x) = X, Q{x) = 1 and Q{x) = 0. 

Remark 6. In equation (12.151) the degenerate case p = const is possible, and in the equation 
(I2.16P the same degeneration leads to a special case of equation (12.141) . 
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Let us prove Theorem 1. We note that the provided below proof contains the algorithm for 
reducing an arbitrary integrable equation (10. 2 p to one of canonical forms (I2.10l) - (l2.23p by point 
transformations. 

Proof. According to Lemma 2, each integrable equation can be written as follows, 
(2.27) ut = u^xx + -42(mx, m, x)ul^ + Ai(m^, m, x)u^.j, + Ao(Mx, m, x). 

It is easy to see that the density of conservation law (12. 4p reads as (12. 7p . where /i = 3 y42,„i— 4 A\. 
Relation (12.81) leads us to two equations 

9^2 , dA2 

9—^ - 36 A2^ + lQAl = 0, 

OUl OUi 

dAo dAo\ „ . rOAo . „ d^Ao „ S^A 



du dx I \ dui ^ ) dxdui dudui 



The first of the equations has a solution of the form 

at that, the second equation becomes 

/ , „ dB dB\ d'^B „ d d'^B 
2.29 2 AiS + 3 — + 3 m— = 3 fi— 

\ ox ou J OUI c'M^; 

In view of the formula for the function yl2, it is clear that without loss of generality the 
leading coefficient of the polynomial B{ui) can be assumed to be one. This is why we have 
three cases, 

1. B = ul + Bi{x,u)ui + Bo{x,u), 11. B = ui + Bo{x,u), III. i? = L 

Equation (I2.29p holds identically in the cases II and III, and in the first case it determines the 
function Ai, 

3 fdB dB' 



2B \ dx du 

Case I. Under the point transformation u = (p{x,v) the function B changes by the rule 

B{x,v) = {ifyVi + (p^f + Bi{x,(f){(fi,vi + (f^) + Bq{x,(p). 

Therefore, taking for Lp any solution to the equation ip^ = —^Bi{x,ip), we reduce the issue to 
the case Bi = 0. 

Returning back to the study of second integrability condition (12.41) . we find that 
d utBo,, ul \d'Ao , 3 f d^Ao d^A 



dV' 4(Mf + 5o)3 6 



Ml^^T^, TTT- + ^iBo,U,Ui 



du\ u\ + Bq \ du\ du\ 

+ Z2U\ + Z1U2 + Zq, 



where the expression $ depends on the derivatives of the functions Bq and vanishes as Bq is 
constant; Zi are some functions in x, u, ui. Equating the coefficient at U2 to zero, we find that 
Bq x = and therefore, B = ul + Bq{u). By an appropriate point transformation u ^{u) we 
convert Bq to a constant Cq being zero or one (case I.l), or zero (case 1.2). 

Equating the coefficient at u\ in (I2.30p . where $ = 0, i?o = cq, to zero, we find the function 
Aq. As a result, equation (I2.27P casts into the form, 

(2.31) Ut = u^^^ - + Aq{ui, u, x), 

2{ui + Co) 
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where Aq is defined by one of tlie following formulas, 

1.1. Co = 1, Aq = ao{ul + 1)^/^ + aiUi{ul + 1) + 02^1 + 03, 

1.2. Co = 0, Aq = — + fliM? + a2Ui + 03. 

Ml 

In both cases = ai{x, u). 

In the case I.l it follows from the further implications of the second integrability condi- 
tion that Co, 02, and 03 are constant and the function ai is independent of u. Moreover, 
a'l = —Saia'i, aoa'i = a^a'^ = 0. If a[ ^ 0, then eliminating the constant 02 by the Galilean 
transformation, we obtain equation f l2.16p . If a[ = 0, then up to the Galilean transformation 
we get equation fl2.14|] . 

In the case 1.2, equating the coefficient at m| in (12.301) to zero, we find the equation 

dai 4 ,da2 3 da2 2 , o'^^s , dao 
5^ "1 ~ 4^ Ui - — Ml + 2— m + — = 0. 
ox ou ox ox ox 

It yields that 02 is constant and the functions ao, ai, and 03 depend on u only. The constant 

02 is eliminated by the Galilean transformation, and one of the functions ao,ai, or 03 can be 

made constant by an appropriate transformation u — ({>{u). 

1.2.1. If ao 7^ 0, then without loss of generality we can assume that ao = 1. In this 
case the second integrability condition is equivalent to three equations, ag = 0, a^a'^ = 0, 
a'C + 8aia[ = 0. If a[ 7^ 0, then letting ai = — 3/2p, we arrive at equation (I2.15p . If a[ = 0, 
then the transformation n — t- n + a^t is admissible and it eliminates the constant a3. In this 
case we get the equation coinciding with (I2.15P for a constant function p. 

1.2.2. If ao = 0, then by a transformation u — ip{u) one can simplify a3 or ai. If = 0, then 
by such a transformation one can eliminate ai, and we obtain equation (I2.17p . If a3 7^ 0, then 
by a transformation u — )■ (p{u) we make a3 constant. Then the second integrability condition 
implies a'l = that allows us to employ the transformation n — )■ n + ast eliminating a3. That 
is, we arrive at the case = considered above. 

In the case I the complete classification has been obtained just by conditions (12. 3p (Lemma 
2) and (12. 4p . It has happened to be possible because pi is the density of high (second) order. 

Case II. In this case B = Ux + Bq{x, u). By a transformation u — ?■ ip{x, u) one can eliminate 
the function Bq. Letting Bq = 0, we find that 

It is easy to check that 

d u\ ( 9Mi 9Mi 



m"^ ~ a + + - ^) + H- 0(2), 

where hi and g are some functions of Mi, m, x. Equating to zero the first terms in these 
expressions, we obtain a system being reduced to one second order equation 



392^1 Ml , 

OU\ OUi 



It yields Ai = ai(x, u)ui + a2(x, u)^/Ul, and equation (12.271) reads as 

3^2 

(2.34) ut = Uxxx - -r^ + (aiUx + a2^u'^)un,x + ^o(mx, m, a;). 

4M.r 
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Equating the term at in fl2.32p to zero, we obtain two relations 

(2.35) 3—— = aia2, 3— h a, = 0- 

ou ox 

For the equations of the form (12.341) we have 

(2.36) ^pi ~ Zsul + Z2ul + Z1U2 + Zo, 

where Zi = Zi{ui,u,x). Equating the expression Z^ to zero, we obtain a hnear inhomogeneous 
fourth order equation for Aq, which determines the dependence of the function Aq on Ui, 

Ao = -Ui 6— h 0-1 ]+-aia2Ui' + a^ui + a^Ui ' +05^1 + 06, 

9 V aw / 3 

where ). Now the dependence of all the coefficients Ai on ui is determined and this 

is why one can split the equations w.r.t. ui under the integrability conditions. For instance, 
the coefficient Z2 at in (I2.36P is linear in ui, and this is why the identity Z2 = leads to 
two identities. These equations read as 

(2.37) 3—— = aia4 + 2a2^:—, 6— 30204-12^:^ h2a2^r— = 0. 

ou ox ou ox ox 

The splitting w.r.t. Ui in conditions (I2.3l) - (l2.5p in view of equations (I2.35P and (I2.37P yields 
several additional equations. The simplest of them are 

(2.38) 0203 = 0, a2 3— 0204 = 0, 



du 

(2.39) 27^-18ai^ + 2a^ = 0, 

ow' ou 

dap. da-i [ dao \ 

2.40 IT^ = 0, 3-^ = 2 aiag, 02 2-^ - aA = Q. 

Ox Ou \ Ox J 

To analyze equations (I2.35p - (l2.40p . it is natural to consider two cases, II. 1 02 = or II. 2 
02 ^ 0. 

II. 1. If a2 = 0, then it follows from (I2.35P that Ci = ai{u). The point transformation 
u —J- (p{u), where (p satisfies the equation 3(p" + 2 (ip'^ai^ip) = 0, vanishes ai. Taking into 
consideration relations (I2.35I) - (I2.40I) . we can write equation (I2.34p as 

(2.41) ut = u^xx - + a3{x)ul + a4{x)ul^^ + a5(x, u)ur, + aei^u), 

4m^ 

where 05 = a{x) + 2ua'^{x). For this equation conditions (I2.3l) - (l2.6p are equivalent to simple 
relations 

03 = C3, 04 = C4, as = a(a;), Uq = ci + C2U, C3C2 = C4C2 = 0, c^a' = c^a' = 0, 

where q are arbitrary constants, a is an arbitrary function. 

If C3 7^ or C4 7^ 0, then 05 and ag are constants which can be eliminated by the transfor- 
mation X ^ X + a^t, u u + a^t. As a result we get equation (12.181) . If 03 = 04 = 0, then the 
equation reads as 

3u^ 

Ut = Uxxx - -r^ + Oi{x)u^ + Ci + C2U. 

If here C2 = 0, then the transformation u ^ u + Cit eliminates the constant Ci. If C2 7^ 0, then 
by the shift u ^ u — Ci/c2 we vanish ci, and then by the transformation u — uexp(c2t) we 
eliminate C2. Thus, in all cases we arrive at equation (I2.19p . 
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II. 2. In view of 02 7^ we set 02 = (3/-\/2) exp(^/^/2), then equations fl2.35p are reduced to 
ai = ?)/2%l)u and to Liouville equation 

(2.42) + = 0. 

Then from (|2:38|) and (ICTj) we get 

3 e^'"^ dip da5 3 e"^ dip 



(2.43) 03 = 0, ag = aQ{u), 



-\/2 dx' du 2 dx^ 



and equation (12.391) is reduced to fl2.42l) . Besides the aforementioned relations, conditions 
f l2.3p - ( 12. 6 p imply exactly one more equation, 



da2 \ ^ „ 9^02 



dx J dx 



0. 



The form of equation fl2.42p slightly differs from the standard one, and this is why we provide 
its solution 

, , 2aV' , , , 

ip = lii- --, a z/ ^ 0. 

[a — vy 

It yields the formula for 02, which we write as 



Ja'{xy{u) 

02 = 6- . 

a — V 

Now it is easy to find ai, 04, and 05, 

3 a" 3a'2 
^5 = 7 + g(x), 

where q is an arbitrary function. Substituting 02 and 05 into equation (12.440 . we find the 
functions and g, 

Co + ciz/ + C2Z/^ 3 a"^ a'" cq + cia + C2a^ 



oe = : , q 



a'2 



where Cj are the constants in the separation of variables. 

Since v' ^ 0, by the point transformation u = z/(u) the result can be somewhat simplified 
and we arrive at equation (12. 20 p . 

Case III. As -B = 1, it follows from (I2.28P that A2 = 0. Then the total derivative of po w.r.t. 
t casts into the form 



Equating the expressions at U2 and to zero, we obtain respectively 

Ai = ao(x, u) + ai(x, u)ui + a2{x, u)ul 

and 

(9a2 da2 2\ 

"9i7 ~^ ~dx ~ ^"^^^^ ~^ '^^^^ ~^ «2^ti) = 0. 
Splitting the last relation w.r.t. Ui, we find that 02 = 0. 
Thus, equation (I2.27P becomes 

(2.45) ut = u^xx + {ao{x,u) + ai{x,u) Ux) Uxx + Aq{ux,u,x). 
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For this equation by a transformation u — )■ (f{x,u) we can reduce the issue to ai = 0. After 
this simphfication one can get easily that 



dt^^ duf ~'~ \ ^M^wf ~'~ dxdu\ Sn^ J ^ 

By this, equating the coefficient at u\ to zero, we find 

Aq = 02^1 + a^ul + 04^1 + as, 

where ). Then the coefficient at gives an equation, which can be sphtted w.r.t. 

Ui into two following ones, 

2.46 17^ = 0' 2aoa2 = 3^. 

Returning back to the analysis of conditions fl2.3p -( p3|) . we can split now also w.r.t. Ui. In 
view of f l2.46p it allows us, in particular, to obtain 

(2 47) = = a2 (3— - 2—^ = 

Consider alternative cases III.l 02 7^ or III. 2 02 = 0. 

III.l. By the transformation u — )■ u fi{x) we can normalize a2{x), a2 = —1/2, which yields 
ao = 0, 03 = a{x). Then by condition fl2.4l) we find 



a4 = fo{x) + cie^- + C2e-2^ as = - -(4 ««' + 3 



6 

where Ci and C2 are constant and, in addition, Qa = q/i = Cj/q = 0, i = 1,2. 

If Ci 7^ or C2 7^ 0, we obtain the equation differing from (12.131) by the Galilean transforma- 
tion. If Ci = C2 = 0, then condition (12. 6p determines /o and /i, 

2 2 4 

fo = ki--(y'^, fi = k2- -{kia + a") + —a^, 

2 /" 

where ki and /c2 are constant. Making the transformation n — n + - / a{x)dx in the obtained 

equation, we reduce it to the form Ut = U3 — ul/2 + kiUi + k2 being equivalent to a partial case 
of equation (I2.13p . 

III. 2. Equations (I2.47P imply oq = bi{x)u^ + b2{x)u + b3{x), = b4{x)u + b^{x). In this case 
the equation (12.451) is simplified by the transformation u — u fi{x) + f2{x)- There are three 
non-equivalent cases III. 2. a ao = 3u^ + b{x), III.2.b ao = 3n and III.2.C ao = 0. 

In the first two cases a simple check of conditions (I2.3I) - (I2.6I) leads us to equations (I2.2ip and 
(I2.22p . respectively. 

In the case III.2.C the form of the equation is determined by three conditions (I2.3p - (l2.5p . 

ut = u^xx + a2,{x)ul + a4(x, u)u^ + a5(a;, u), 

where a4 = biu^ + 62 w + 63, as = b^u^ + b^v^ + b^u + 67, bi = bi{x), 1 ^ i ^ 7. Integrability 
conditions (I2.3p - (l2.6p yield a cumbersome system of equations for the functions bi, studying of 
which yields several forks. 

1. If a3 7^ 0, then a3 = co, 61 = 62 = &4 = &5 = ^6 = 0, 67 = cix + C2 + ^b'^ + ^63. The 

transformation u ^ u / b^dx yields the equation ut = + cou^ + cix + C2. The sixth 



2 „ 

integrability condition implies Ci = 0; then by the transformation m — )■ m + C2t we eliminate C2 
and obtain flZT^ . 

2. If a3 = 0, then bi = const, 64 = 0. Then again forks appear; 
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2.1. If bi 7^ 0, then 63, 65, Bq, and 67 are expressed in terms of 62(2;) so that the transformation 
n — 7- n — 62/ (2&1) leads to the equation equivalent to (12. lip . 

2.2. If bi = 0, we obtain 62 = const, 62(^6 - ^s) = 0, 62(^3" + ^3^3 - ^2^7) = 0. If 62 7^ 0, 
the transformation m — )■ u — 63/62 leads us to the equation equivalent to (12.101) . Otherwise we 
obtain linear equation (12.231) . 

2.3. Comments to the list of integrable equations. It was shown in the previous section 
that each integrable equation (10. 2 p is reduced to one of equations (I2.10l) - (l2.23p by a chain of 
point transformations. Although the answer as the list is not invariant w.r.t. point transfor- 
mations, integrability conditions (12. ip . (12. 2p are so. This is why to check the integrability of a 
given equation it is not necessary to reduce the equation to one in the list. According to the 
proof of Theorem 1, it is sufficient to check four conditions (12. 30 - 02. 61) if the equation belongs 
to the classes I or II, and six conditions (12. 2p if the equation belongs to the class III. It can 
be shown that if the right hand side of equation (10. 2p is independent explicitly of x, then for 
the equations of the class III it is sufficient to check conditions (I2.3p -( l2l6|) . 

The orders of canonical conservation laws (I2.ip . (12. 2 p are discrete invariants of the group 
of point transformations. The analysis of the structure of these conservation laws shows that 
these equations are divided into two groups. For the first group of the equations (let us call 
them S'-integrabl^) the canonical series contains conservation laws of arbitrarily high order. 
We observe that this property is more restrictive than just the condition of the existence of an 
infinite series of conservation laws for the equation. For instance, the linear equation ut = Uxxx 
possesses an infinite series of the conservation laws with the densities u\, /c G N. However, all 
its canonical laws are trivial. 

For the equations of the second group (C-integrable equations) among the canonical conser- 
vation laws there are only a few non-trivial ones. Equations (I2.19p - (l2.23p are S-integrable. The 
equations (I2.19P and (I2.23P have no non-trivial canonical conservation laws. Equation (I2.20p 

has just one nontrivial canonical conservation law of the first order po ~ ^ ^ . Each 

u — a[x) 

of equations (I2.2ip and (I2.22p has one canonical conservation law of zero order, po ~ and 
Po u, respectively. 

All canonical conservation laws with even indexes of S-integrable equations are trivial [7], 
and the orders of odd laws increase by one, but the initial orders in these sequences are different. 
In Table 1 we provide the order of first four odd canonical conservation laws for all S-integrable 
equations. 

For equation (12.180 we provide the orders of the densities in the case of generic constants. 
In the case Ci = we have pi ~ 0, and other orders remain to be unchanged. If C2 = 0, the 
orders are equal to 1, (p3 ~ 0), 2, and 3. 

Remark 7. Equations (I2.10p - (l2.18p are integrable by the method of the inverse scattering 
problem while (I2.19p -f l2.22p are linearizable by the differential substitutions (see Section 2.4). 

Table 1. Orders of canonical conservation laws. For zero order conservation laws we indicate in the 
brackets to what the density is equivalent 



pi 


fl2.10p 


fl2.1ip 


fl2.12p 


02.131) 


02.14P 


02.16P 


02.15P 


02.17P 


02.18P 


pi 


0, (~ u) 


0, (~ u') 


0, (~ 0) 


1 


2 


2 


2 


2 


1 


P3 


0, (~ u') 


1 


1 


2 


3 


3 


3 


3 


2 


P5 


1 


2 


2 


3 


4 


4 


4 


4 


3 


P7 


2 


3 


3 


4 


5 


5 


5 


5 


4 



The terminology is due to F. Calogero. 
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If in the formulation of the original classification problem one assumes that the right hand side 
of equation (10.21) is independent explicitly of x, then the answer changes just for C-integrable 
equations. For equations f l2.19p . fl2.2ip . f l2.23p arbitrary functions are replaced by arbitrary 
constants and then these constants can be eliminated by point transformations. 

Formula fl2.2Qp contains two C-integrable equations independent of x explicitly, 

3 ^^2 

Ut = Uxxx 7 — 'iUxxU~^iVUx + 1 + Ux + 1) 

(2.48) Aux + 1 

+ 6u-\x{ux + 1)^/^ + 3u-\x {Ux + l){Ux + 2), 



_ _ 3 _ Uxx jUx + l) COshw UxxVUx + 1 

Ut — Uxxx , , -, <J .1 +0 



(2.49) 



4 + 1 sinh u sinh u 

Ui.(Ui + 1)3/2 cosh U Ux{Ux^\){Ux^1) 2/ , 

- b T-Ts + 3 —2 + + 3). 

smh u smh u 

Equation fl2.48p can be obtained from (I2.20p with a = x, Cq = Ci = C2 = by the transfor- 
mation n — 7- n + X. Equation (12.491) is obtained in the case a = e^^, cq = ci = C2 = by the 
transformation u — > e'^^'^~^^\ It is impossible to obtain from (I2.20p any other equation indepen- 
dent on X. It follows, for instance, from the results of independent classification of equations 
(10. 2 p without explicit dependence on x. 

The situation with equation (I2.22p is rather instructive. Let the functions 7 and /3 be 
constant. Then the constant 7 is eliminated by the Galilean transformation. It is easy to check 
that if /3 7^ 0, then the integrability conditions hold true but the canonical conservation laws 
depend explicitly on x. It is impossible if in the classification of equations possessing generalized 
symmetries we restrict the generalized symmetries not to depend explicitly on x. In this case 
/3 = and equation (12.220 is just a third order symmetry for the Burgers equation. If the 
dependence on x for the symmetries is admitted, then the constant /3 in the answer should be 
preserved. 

2.4. Differential substitutions relating the equations in the list. We say that a differ- 
ential substitution 

(2.50) M = $(x, M, Ml, . . . , Mfc) 

acts from the equation 

(2.51) Ut=Un + g{x,u,Ux,...,Un-i) 
into the equation 

(2.52) Ut = Un + f{x, U, Ux, • • • , Un-l) 

if for each solution u{x,t) of equation (I2.5ip formula (I2.50p provides a solution for equation 
(I2.52p . The number k is called an order of the substitution. Since for k > transformation 
(I2.50p has no inverse transformation of the same form, equations (I2.5ip and (I2.52p in this 
definition are not equal in rights. If a differential substitution reads as (I2.50p . where u satisfies 
an equation (A) and u does equation (B), we express it as (B) — (A). 

The most known differential substitution is the Miura transformation u = Ux — u^ relating 
the Korteweg-de Vries equation 

Ut = Uxxx + QuUx 
and the modified Korteweg-de Vries equation 

Ut = Uxxx -QU^ Ux- 
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Other substitutions relating the main equations of the hst were found in [52]. The issue on the 
invertibihty of the differential substitutions was treated in [53] . 

The orders of possible differential substitutions relating S'-integrable equations can be found 
by Table 1. Namely, if equations f l2.5ip and f l2.52p are related by the substitution fl2.50p . then 
the orders of the canonical conservation laws with sufficiently large indices for fl2.5ip are greater 
by k than the orders of the canonical conservation laws with the same indices for fl2.52p . For 
instance, if equations fl2.16p and fl2.10p are related by a differential substitution, then it acts 
from (12.161) into fl2.10p and is of the third order. 

Below we provide the information on differential substitutions relating various integrable 
equations in the list. Since the superposition of differential substitution fl2.50p and the Galilean 
transformation leads out of the class of the substitutions of the form fl2.50p . sometimes in order 
to find a substitution it is necessary to add a term cUx to the right hand side of equation fl2.5ip . 
All such cases are mentioned below in the text. 

I. S'-integrable equations. It happens that for all S'-integrable equations except the 
Krichever-Novikov equation fl2.15p there exist substitutions acting in the Korteweg-de Vries 
equation. For equation (12.151) such substitution exists only when the Weierstrass function 
degenerates or, which is the same, the polynomial Q in formula (12.240 has multiple roots. 

Let us provide all substitutions in the Korteweg-de Vries equation. If from a given equation 
there exists several substitutions into (I2.10p . we provide all of them. We replace equations 
(I2.15P and (12.161) by (I2.24p and (I2.25p . respectively, since after that the substitutions look more 
compact. 

( 12.110 — 7- ( 12. lOh : u = ±i\/6ui + u"^ + X; at that Ut = u^xx + u^Ux + Xu^- 
( 12J^ ^( I2JX)D : u = 2ui. 

( BTTBD ^d^TTUD: u = 3u2-^ul + 2V^6^ui e"" + cie^" + 026"^". 

(EUlD^dSaOD: u = ^_ - , , , - ^ == + 6con2 + 3aini+ 




3 ai Ml v/ + 1 , where cq = ^/ (cti — a.2)/2 . 



u=^j6 ——^ -^=== j - g^2(g + ,2) + 2^ > where 

Q(a) = 0. 

■"2 , ^ 3 



(EUSD^dMOD: S = v/^3^— ^ + 2c2Ui + -Civ/^. 

Jui 2 



The above substitutions of higher orders are the superpositions of first order substitutions. 
These substitutions relate some of S-integrable equations. The first order substitutions are 
drawn on the graph (Fig. 1). 

The arrows of the graph correspond to the following substitutions, 

( BTT^ ^ (1^7151) : M = In (ui + y^l + y At that, in the equation (Oip there should be 
3 

an additional term -02^1. The constants in the equations are related by the formulas ci = 

3 3 

-(fli + 02), C2 = -(02 - ai)- 
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(I235D^ (12331): u = \n + 



ln(ao + 2 fli n + 02 u^). At that, the polynomial 



Q is written in the factorized form, Q{u) 



2kiu 



k2U^). More- 



(ao + 2ain + a2 u^){ko 

over, in equation fl2.25l) there should be an additional term -(09^2 + ^2^0 ~ 2 aiA;i)ui, and 

the constants Ci and C2 in equation (12.131) are determined by the formulas Ci = -(0002 — a\), 
3 

C2 = -(^0^2 - kl). 



(I233D^(I23TD: n = ±-V6ni + v^e" + 
-y/c2e~", at that, in equation (12.131) there 
should be an additional term 1^C\C2 U\. 
dSHSD^dSmi): n = a + 6^^, h ^ 0, 

where Ci = - a6, C2 = -l? . At that, in 

3 2 
equation (I2.18P there should be an addi- 
tional term c?u\. 

If in equation (12.241) the polynomial Q 
has multiple roots, then there exists the fol- 
lowing substitutions from this equation not 
indicated in the graph, 

U2 ^ Vlux\ 3^2 Q 



1) 

Q 



u 



d 

dx \ Ui u 

{u - of (co + ClU + C2U^) . 



3n2 
2^ 



u 



2 ' 




Figure 1. Graph of the substitutions for 
S-integrable third order equations 



2) u 



d f ^'^2 
dx \ Ui 



I2h\ ?,ul 



— h = co + ciu + C2U^, Q = Qh'^ 



u 



-C2/i^ ci = -(4 0002 - al). 



ui J 2u\ 

(EISlD^dSIISD: = Inwi - In/i, h = ao + aiu + a2U^, Q 

In the case Q = equation (I2.24p coincides with Schwarz-Korteweg-de Vries equation (I2.17p . 
Equation (I2.17P is related with the Korteweg-de Vries equations by three different substitutions, 

o . " ' ]; 3) S = -3 — + ■ 



(IlalD^dSaOl): l)n = 3^ 

Ui 



^; 2) S = -3 — , 
2ul' ' ui 2u{ 



3 , 3 MS I U2 U\ 



ui 2 u{ 



u 



All of them are superpositions of first order substitutions. Besides aforementioned, in these 
superpositions the following substitutions from equation (I2.17P are involved, 
dSIUD^dSIIS]): l)u = In(Mi), ci = C2 = 0, 2) u = In(Mi) - \n{u^ + Ci/6), C2 = 0. One more 
substitution is obtained from 2) by the replacement u — ?■ — u, C2 ^ ci. 
II. C-integrable equations. 

( 12J^ ^ ( 127231) : M = 017, at that, in equation fl2:23|) /3 = 
dSOlD^dMlD: n = v/^i/(2n). 

(127231) ^ (12:221) : u = Ui/u, at that, in equation fl2:23|) a 



2 

: 7. 
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(I2.20p — 7-( l2.22p : u = -{a' + y/uxa' ) — - a"{a') ^, where ^ = a{x) — u. At that, in equation 

1 3 

7 = - - ^ (a'7«')^ - (co + Cia + C2a^)/a\ 

13 = ^ a(^) («')"' -2a"a'"{a')-^ + ^ {a" /a'f - ^ a\a')-\cQ + cia + Caa^). 

3. Fifth order equations 

In this section we find all equations fl0.4p possessing infinite sequences of local generalized 
symmetries. In the process of classification we employed necessary integrability conditions 
which follow from the existence of formal symmetry [HI?] and are written as canonical conser- 
vation laws. It is easy to check that each third order integrable equation (I2.10p -( l2.23p has a 
fifth order symmetry fl0.4l) . It turns out that if one exclude these symmetries from the consid- 
eration, then the list of other integrable equations coincides (up to the equivalency) with the 
list obtained in works [9l[T8] , where the equations possessing generahzed conservation laws were 
studied. 

In Subsection 3.1 we provide the complete list of integrable equations flU.4l) not being sym- 
metries of lower order equations. The equations in the list differ slightly by the form from 
equivalent equations in [9l|18]. Subsection 3.2 contains a new recurrent formula for the inte- 
grability conditions. We note that in the works [HllIS] only several simplest conditions were 
provided. For a given equation (10.41) the integrability conditions can be easily checked one by 
one by computer. 

In Subsection 3.3 we adduce a schematic proof of the classification theorem. It contains 
an algorithm for reducing an integrable equation (10. 4p to one of the canonical forms from 
Subsection 3.1 by point transformations (ll.24l) - (ll.26l) . In other words, once we employ point 
transformations, we indicate what we normalize by them. We hope that following this algorithm 
and guidelines given in the text, the interested reader can easily recover all the details of rather 
hard calculations. Post factum it follows from the proof that if equation (10.41) satisfies first 
ten integrability conditions, then it is integrable. We observe that the equation ut = + uui 
satisfies fist nine integrability condition, but does not tenth. 

3.1. List of integrable equations. 

Theorem 2. Suppose nonlinear equation (10.41) satisfies two conditions, 
1) there exists an infinite sequence of generalized symmetries 

(3.1) Ur^ = Gi{u, ...,Un,), Z = l,2,..., Ui+i > Hi > ■ ■ ■ > 5] 



2) there exist no symmetries (13. ip of orders 1 < < 5. 
Then the equation is equivalent to one in the list 



(3.2) 


Ut 


= M5 - 


h 5uus 4- 


- 5uiU2 + bu^ui, 


(3.3) 


Ut 


= Ms - 


h 5uus 4- 


25 , 2 

- —U1U2 + 5m Ml, 


(3.4) 


Ut 


= Us - 


h 5nin3 - 


5 3 


(3.5) 


Ut 


= Ms - 


h 5uiUs - 


^ T ' ^ 3""'' 


(3.6) 


Ut 


= Us - 


h5(ni - 


u'^)u3 + 5^2 — 2OUU1U2 — 5ul + 5u'^Ui 
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= M5 + 5(m2 — Mi)m3 — 5uiul + ul, 
Ut = U5 + 5{u2 -ul + Aie^" - A^e"^") - Sum^ + 15(Aie2" + 4A^e-^") 

+ - 90A^e-^" ul + 5(Aie2" - A^e"^")' 
Ut = U5 + 5{u2 - Ml - A^e^" + A2e"'") % - 5miM2 - 15A^e^"MiM2 

+ u? + 5(A?e2"- A2e-")2ni, A2 7^ 0, 

M2M4 , ^ujus , p. //il , 2^ c: /^/^l , \ 2 

Ut = u^- 5 h 5 — 5- + 5 h yU2n, U3 - 5 ^ + ;U2ni 

Ui U{ \Ui J \u{ J 

^2 

- 5— + 5/iiyU2Ui + yUo^i, 

U2U4 15 u\ 65 ^2^3 / /ii 2\ 135^2 

Mj = U5 - 5 h 2~ + ^ ^ /"S^il ^^3 - -77 3 

ui 4 Ui 4 nf \ui J 16 nf 



/7/ii /i2Mi\ 2 , p. 2,25 

\4m| 2 y Ui 



5 U2U4 5 ul r-^lus 5U2U3 1/2 , 2\ 35^2 

2 Ml 4 Ml 2y7ll 16 UI 



3 ' "V4wi ^AI 47 ^ ' 3 ^ ^ ' ' ^ ' 3 



5 / — Ml 5 2/— Ml 2 ^ / r\2 
= Ms + 2 j2 "2M4 + ^ M3 + 5 /i (Ml + /) M3 

5 4mi2 -8M1/ + /2 2 5 2-9^3 + 18m?/ . 
+ 4 ^ "^"^+16 P 

, 5/i (4/-3mi)(mi + /)^ 2,2/ , n^foff , i^ 
+ ^ Y2 U2 + {ui + f) [2f{ui + f) -1), 



5 / — Ui 5 2/ — Iti 2 / r2 2\ 

Mt = M5 + - ^2 ^2^4 + - ^3 - 5 W (/ + n JU3 



54M?-8ni/ + /2 2 5 2-9m? + 18m2/ 
+ i J. -2-3 + ^ ^ M2 

5 5n;-2M^/-llMi/^- 2 2 5 2 . , ^ , ,2^ 

+ - W — U2- -UJ (Ml - 2 m/ + 5 / )MiM2 

+ 5a;V/'(3Mi + /)(/-Mi), 

, 5 / -Ml , 5 2/-M1 2 ^ 5 4m?-8mi/ + /2 2 

Mi = M5+ 2— ^«2M4+^— ^% + - U2U3 

5 2-9mi3 + 18mi2/ , 2mi3 + mi2/-2mi/2 + 1 2 

+ 16 p U2 + 5UJ M2 

- 10 a; M3(3 Ml/ + 2 M? + 2 f) - 10 to' {2 f + mJ + mi^) M1M2 
+ 20a;VK^- l)C«i + 2/), 
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+ 4 ^ ^2-3 + ^ 



(3_1Q) -I0ui3uj + 2u, + 2f )u3--c 

2ul + ulf-2mp + l 2 , ul + 5p-2uj 

+ 5UJ Uo + OCU 5 U1U2 

- 10 {2f + uif + Ui^)uiU2 + 2{]uj\i{u\-l){ui + 2f) 

Here Ai, A2, /i, Hi, /i2, and c are parameters, the function f{ui) solves the algebraic equation 

(3.17) (/ + Mi)2(2/-Mi) + 1 = 0, 
and uj{u) is any nonconstant solution to the differential equation 

(3.18) uj''^ = Auj^ + c. 

Remark 8. All the equations in the list of Theorem 2 are S'-integrable. In the proof of the 
theorem it has been established that each C-integrable equation (10.41) is a symmetry of some 
third order C-integrable equation in list fl2.10p - fl2.23p . 

Remark 9. If Ai = A2 = in equation (13.81) . it coincides with (13. 7p . If Ai = A2 = in 
equation (13.91) . it also coincides with (13. 7p . and if A2 = 0, then (13.91) coincides with (13. 8 p as 
A2 = up to the change Ai — t- — A^ in the latter. 

Remark 10. If /i2 7^ in equation (I3.10p . the substitution u = \nv reduces it to 

(3.19) Vt = V5- 5 h 5 — 5- + 5yU 5- - — , 

vi v{ \ v\ v{ Vx 



where [i = [i\C, c = ^—[i2- 

Remark 11. If /i2 7^ in equation (13. lip , the transformation u = c~^lnw reduces it to 

3.20 vt = V5- 5^ - -r^ + - + 5/i ^ + -V2 - - 5/i'-, 



where /x = /iic, c = 2y/—fi2. 

3.2. Integrability conditions. The following conditions can be borrowed from the works 

peg. 

Lemma 3. For nonlinear integrable equations ( [0.-^[ ) the first four integrability conditions can 
be written as 

(3 21) '^^ 

dt dui dx 

(3 22) ^ (2(^^y 

dt I \du4^/ du-i j dx^^' 

d [ dF dF dF . f dFV\ d 



^'^''^'^^ dt \^^du-idu4^ du2 ^ydui) j dx^"^' 
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(3.24) 



d 
It 



25 



d dF 
dx duA 



dF 



5—— + 10— _ 



dud V dx du 



dF 

duA 



25 



du2, du4 

2 



+ 



dF 
du?. 



125 



dF 

dui 



dx 



where — is the operator of the total derivative w.r.t. x, and — is the evolution derivation 

dx dt 

w.r.t. t by virtue of equation 



The given conditions follow from the existence of a formal symmetry. But technically it is 
more convenient to employ the described in Appendix 3 method of calculating the densities of 
canonical conservation laws 



(3.25) 



d d 

(jjV (JjJb 



by the logarithmic derivative of a formal eigenfunction of the linearization operator for equation 
(10. 4p . This approach allows us to obtain the recurrent formula 

Pn+4 =-(^n — - [F'uo^nfi + -^ui<^n,-l + F^^^n 2 + -^«3^n,-3 + Fu4^n-A + FuiPn\ 

5 5 



(3.26) 



n+3 71+2 ^ n n+1 / ^ \ 

2 PiPj - 2 XI APjPfe - 5 P^PjPkPlPrn + Yl ( d^P' ) 

^ ^ 



d 

dx 



Pi 



d 

dx 



Pn + '^Pn+l + X PiPj 



n+1 



(p d 3 d 

Pn + 3— p„+l + 3p„+2 + 1^-1-Y P'P^ + ^ X] P'P^ + 5Z P'P'P'' 



dx 



^F 



U4 



n+1 



d^ (f d d^ d 

IZ^P- + 4T;:^Pn+i + 6;7-Pn+2 + 4p„+3 + 2;73 E PiP^ + E + 



dx"^ 

n+1 



n+2 n+1 ^ f d \ d d ^ " 

+6 J2 PiPj + 4 J] p^PiPfc - Yl ( ) + ^rf;^ S + pipjpkpi 

^ ^ 



d_ 

dx 



1 





c/2 



n+1 



n+2 



rPri 



5 dx^ dx"^ 
2 c/ 



d d"^ d 

00 

n+1 



2 c? X ^ 1 X ^ f d \ d \ 1 v ^ 

+2p„+3 + 3^2^ PiPjPfc + 2 V rf;^^' j di^^ + 2 ^^^^'^^ ^ 2 ^ PiPjPkpi 
^ ^ . 

where n ^ —4, = 0, V/c < 0, is the Kronecker delta, F^. = dF/dui. In formula (I3.26P the 
notation 

n 

Oifej . . . = X aihj... pz 



h^=n 

i^m, j^m,..., z^m 



was employed for multiple sums. The summation indices in formula (I3.26P are non-negative. 

Recurrent formula fl3.26p is published at the first time. It is easy to check that the first four 
integrability conditions in sequence (I3.25P are equivalent to conditions (I3.2ip - (I3.24p . 
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Conditions fl3.26p . ( I3.25P can be used more efficiently for the classification if one ffist studies 
the structure of the densities of local conservation laws for equations (10 ■4p . We recall that the 
symbol 0{n) indicates a function of the differential order at most n. Moreover, we employ the 
symbol Pn{uk) to denote a polynomial of the degree n of the variable Uk, whose coefficients 

have differential order less than k. In what follows we employ regularly the equivalence f—g ~ 

d d 
—g—f implied by —ifg) ~ 0. For instance, we have 
dx dx 

i=0 * 

which yields, in particular, Un+iO{n) ~ 0{n). 

Lemma 4. // a density p of a local conservation law for equation 0^ has differential order 
n ^ 3, then the equation 

(3.27) d_^^2d^dF_ 

dx dul^ 5 (9^2 dui 

holds true. 

Proof. By definition we have 

d dp ( d^ ^ 



dt ^ duk \ dx'^ 

Using the equivalence, we can lower the order of this expression up to n + 2. First we show 
that 

a„ / 



I — n \ / 



k=0 

For this it is sufficient to convert the highest order term. Assuming n ^ 3, we get 
dp d dp d'^p ^, , 

-M„+3 ~ -Mn+23- 7^ = -Mn+2Mn.+17^ h Un+20[n) 



dUn-2 dx dUn^2 dUndUn- 

^ 2 ""-^^Txd^I^d^^ - ""^^rf^^^"^ = + 

Thus, 

d dp dp dp d^ dp d^~^ 

^ ^ dt^ dun ""^^ ~'~ dun-i ^^'^ ~^ dun dx^ ~'~ dun-i dx'^~^ + ( + ) 

Let us convert the ffist term, 



dcf 



dp d^ dp d ^ d'^p 



dUn dx'' dUn dx ^-^ dUndUi 

1=0 



OUndUi dUndUidUj 



1 2 d d'^p d I d'^p d^p 

2 ^'^^dxdul ^'^"^ dx \^ dundui ^'^'^ ^ dunduiduj ^^'^ 

2 fl d d'p d^p ^ d'p \ 
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52 „ n2 . 



The second term in fl3.28p is converted in the same way, 

dcf dp / "^V , d^P 

a2 = Un+4 ~ Un+2 > ^ ^ 5~^i+2 + > , ^ ^i+l^j+l 

- -<+2 + «n+20(n + l)~— 4^<+2 + 0(n + l). 



dUndUn-1 dUndUn-l 

The above calculations are correct if n ^ 2. To convert the remaining two terms in fl3.28p . it 
is important that n ^ 3, 



dcf dp d"- ^ f £ dp\ d 
«3 = T— -7— ^- -7-?!^ -7-;-? ^ 



dun dx^ y^dx"^ dun) dx^~'^ 

d^p ^. fdF 



= ^ ^ + ""-^^^ + ^) + + ^) ^ ^ ^ + + 

a, ^ ^ ^ _ M i:::! ^ ^ _ m (dF_ ^^^^ ^ ^ 

dun-i dx^"^ \dxdun~i) \dxdun~ij \dui " 

= Un+20(n + 1) + 0(n + 1) ~ 0{n + 1). 

Summing up the obtained expressions for ai,. . . , 04, we find 

d 2 (d'p dF 5 d d'p\ 

dtP^""-^' [Md^r^Txd^J^^^''^^^- 

Since a quadratic in higher derivative expression can not be a total derivative of any function, 
we obtain fl3.27l) . 

Corollary. // we have n > 3 in Lemma 4, then the density p is quadratic in Un- Indeed, 
the left hand side of equation 1^3.21^ contains the term PunUnu„Un+i whose differential order is 
greater than 4 if n > 2>. The order of other terms is less and hence Pu„u„u„ = 0. 

We apply the obtained result to the classification of equations fl0.4l) . 

Lemma 5. Assume equation satisfies condition ^3.21\) . Then the function F is qua- 

dratic in U4. 

Proof. Applying Corollary of Lemma 4 to the canonical density p = F^^, we obtain 

dF 

= fi + hui + f3ul, 

where the functions /i , /2 and /s are independent of U4. Substituting this expression into 
( KTHf . we find 

d 2 

^/s = ^ fsifi + /2M4 + fsuD- 

The left hand side of this equation is linear in u^, and the right hand side is quadratic; thus, 
/s = 0. It yields F = fo + fiu^ + - f2u1, where the functions fi are independent of U4. 
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3.3. Scheme of proof of main theorem. 

Lemma 6. Suppose equation Ili0.4\ ) satisfies integrability conditions \3.21\) . \3. 22\) . and 
lis. 23\) . Then the function F is linear in u^^. 

Proof. According to Lemma 5, the function F is quadratic in ti4, F = fo + fiu^ + - f2u1- It 
can be easily verified that 

P2 ~ ul f2 (16 fl - 15 + Z^ul + 0(3). 



In accordance with Corollary of Lemma 4, a cubic in M4 term should vanish and therefore 

fl 



df2 16 ,2 



dus 15 
In view of this equation we find 

Pi^ulf^ + 0{3). 
For this density relation (I3.27P read as follows, 

^/2 = ^/2(/l+/2«4). 



It implies the equation 



du, "5^2' 



which together with the previous one yields /2 = 0. 

Thus, if integrability conditions f l3.2ip - fl3.23p hold true, then equation (10. 4p reads as 

(3.29) ut = u^ + U4fi{u, Ul, M2, U3) + /o(m, Ul, U2, U-s). 

Lemma 7. // a function of third differential order p{u, Ui, U2,u^) is a density of a conserva- 
tion law for equation 113. 29\) . then it is at most quadratic in U3. 
Proof. Letting n = 3 and F = fo + fiu^ in (I3.27P n = 3, we obtain 

(3.30) ^9V^2 9V^^_ 

dx du\ 5 du\ 

Taking into account that fi is independent of u^, we find that Puauaua = 0. 

Corollary 1. The function fi in 113. 29\) is linear in u^. 
Proof. Indeed, it follows from F = /0 + /1M4 and fl3.2ip that fi is the density of a conservation 
law for equation (I3.29p . Thus, as it was proven above, this function reads as fi = gi+g2U3+g3ul, 
where gi = gi{u,ui,U2). Substituting this expression into (I3.30p instead of p, we obtain g^ = 0. 

Corollary 2. If fi = gi + 5'2^3; Oi = gi{u,Ui,U2) in equation h3. 29\) and this equation 
possesses a conservation law with a density p of second differential order, then the equation 

d d'^p 2 d'^p 

holds true. 

Proof. The statement can be easily checked by straightforward calculations. 

Proposition 1. If integrability conditions ^3.21\) - f3.23\) hold true, then equation ( [0.^[ ) reads 

as 

(3.32) Ut = U5 + AiU2U4 + A2U4 + A3ul + {A4ul + A5U2 + AG)u3 + A7U2 + A8ul + Agul + AioU2 + Aii, 

where Ai = Ai(u,ui). 
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Proof. By Corollary 1 of Lemma 7, equation (10. 4p reads as fl3.29p . where fi = gi{u,ui,U2) + 
g2{u,ui,U2)u3. Consider then a density pi of conservation law f l3.22p . It is easy to check that 

2 ,2 ^ dh ^ dh dfo 

Pi ~ ^ A + T^Ml + -^—^2 + -W—U^ - -W—- 
5 OUq OUi OU2 OU3 

In this expression all terms except the last one are at most quadratic in u^. By Lemma 7 the 
considered density must be quadratic in M3 and therefore, the function /q is cubic in M3, 

fo = 94 + 95^3 + 96^1 + gjul, Qi = gi{u, ui, U2). 

Due to the obtained results, densities of conservation laws (13.221) and (I3.23P are equivalent 
to the expressions 

,2 f::d92 , . 2 



Pi ~ ui (^5^ + 2 - 15^77 J + 0(2), 

P2 ~ ^4 (5o|^ - 25^ + 30 (72!^ + 89I- 909297^ + P2(%). 

According to Lemma 7, the coefficient at in the second formula should vanish. Moreover, 

the condition — pi ~ leads to extra four equations relating the functions g2, 97 and their 
at 

derivatives w.r.t. U2. By these equations it is easy to obtain g2 = gj = 0. 

Thus, F = giU4 + (74 + g^u^ + geu"^. Now the density of conservation law (I3.2ip equals gi and 
we can substitute p = gi and g2 = into (I3.3ip . 

dd'gi 2 9V 

91- 



dx dul 5 dul 



As above, by this we obtain the linear in higher derivative function gi = Ai{u, Ui)u2 + A2{u, Ui). 
In view of the obtain results, condition (I3.23P yields g^ = A-^^u, ui). Then by condition (I3.22p 

we get = 0. And finally, bearing in mind all the obtained results, we find by condition 



(K23\\ that — ^ = 



9^94 

In studying equation (I3.32p the following lemma will be useful. 

Lemma 8. Equation i\3. 32^) preserves its form under point transformations u = (f{v). 
Some of the formulas for the conversion of the coefficients Ai look simple, 

A,{v) = if'A.iu), A2iv) = A2{u) + ip"vlA,{u) + 5ip"{ipT'vu 

(3.33) As{v) = if'Asiu), A^{v) = ip'''A,{u), Aj{v) = ip'^A-r{u). 

Other formulas are much more bulky and we omit them. 

It can be checked that the first six densities of the canonical conservation laws for equation 
(I3.32P are equivalent to 

(3.34) po = -hAiU2 + A2), pi ~ i?i = z/'iMa + ^2^2 + ^^3, 

5 

(3.35) p2 ~ i?2 = i'iul + tp^uj + i)6U2 + ^7, P3 ~ i?3 = ^8^*3 + ^9^2 + ^10^*2 + • • • , 

(3.36) p4 ~ = 1puU2ul + 1pl2ul + ^13^2 + . . . , 

(3.37) P5 ~ -R5 = i^iul + ipi5ul + {ipmul + tpnU2 + i^is)ul + tpwul + 
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where the coefficients ipk are expressed in terms of the functions Ai and their derivatives. For 
instance, ipi in f l3.34p and f l3.37p reads as 

ipi = — {2Al- 5^4 + 10 



25 V dui 

Lemma 9. // equation l{3. 3^) has a conservation law with a density p{u,Ui, . . . ,Un) of 
differential order n ^ 2, then 

(3.38) p ~ ai{u, ui)ul + a2{u, ui, . . . , Un-i), 

at that, 

dai , dai 

OUi OUq 

For n = 3 and n = 4 the structure of the densities can be easily specified. If n = 3, then 

and ai{Ai - 2 A3) = 0. If n = 4, then 

p ~ aiul + ^2^3 + {ce^ul + '^4^2 + ^5)^3 + (^iu, Ui, U2), ai = ai{u, Ui), 

where /3 is a polynomial of sixth degree in M2. 

Corollary. The coefficients ip^ in 113. 35\} . ipu and ipi^ in 113. 36\} are zero. 

The form of equation (13.321) depends essentially on the orders of its canonical conservation 
laws. Among integrable equations (13.321) . there can be equations of the two following types, 

I. Equations possessing no generalized canonical conservation laws. In other words, all canon- 
ical densities for the equations of the ffist type are equivalent to densities of zero or ffist differ- 
ential order. 

II. Equations possessing generalized canonical conservation laws of orders ^ 2. 

In the case I one should equate all nontrivial terms of higher order in the densities of the 
canonical conservation laws to zero. This is why in expressions (I3.34I) - (I3.37I) there should be 
V'l = '04 = V's = V's = V'g = '010 = ■ ■ • = 0. In particular, 

2 dAs 
= -Al + 2-^. 

5 oui 

By the equation 1^4 = one can express Aj in terms of Ai and ^3, and hj tp^ = one can 
express As in terms of Ai, A2, A3 and A5. From six conditions pi ~ hi{u, ui), i = 1, . . . ,6 one 
can extract a cumbersome system of differential equations for the remaining functions Ai. In 
this system there is the following closed subsystem of the equations for Ai and ^3, 

(3.40) A,.iA., '^y-Al 

The latter of these equations has two solutions Ai = and Ai = — |(mi + a(-u))~^. If a{u) 7^ 0, 
then by point transformation u — )■ ip{u) one can normalize a = 1. Thus, there are three possible 
cases, 

I.a. Ai = 0; I.b. A, = -^u^^; I.e. A, = -^{u, + iy\ 

Case I.a. It follows from the equations ipi = that A2 = gi{u) + g2{u)ui. Employing 
the point transformation u — )■ ^{u), we can assume that (72 = (see (13.331) ). After that all 
remaining functions Ai{u,ui) are happened to be polynomials with constant coefficients. To 
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determine these coefficients we check 10 integrabihty conditions (13.251) . We find out that there 
exist only three integrable equations of the considered type, 

Ut = U5 + U4C1 + C2U3 + C3U2 + C4U1 + C5U + C6, 

Ut = + ^u^Ui + lOn u^{u^ + 4ni) + 25n + 10n2(5n^ + 12n^ni + n^) 

Ut = + 5n 1x4 + lOn u-^ + 15^1^3 + 10^2 + lOn U2 + 50m U1U2 + 5n Ui + 30n + 15n^. 

The second of these equations is a symmetry for equation (I2.21|) . where a = 0. The third 
equation is the symmetry of the Burgers equation Ut = U2 + 2uui (as well as a symmetry of 
equation (12.221) . where /3 = 7 = 0). 

Case I.b. There exists only one integrable equation in this class, 

5U2U4 ^ 5^3 35Ug 



Ut = U5 



5n^ 

-r^ - T7r4 + ku. 



2ui uf Aui IQuf 

It is a symmetry of equation (12.191) as a{x) = c. 

Case I.e. By the integrabihty conditions one can find A2 = f{u){ui + 1) + g{u)^/ui + 1, 
where / and g are arbitrary functions. If (7 = 0, all the functions Ai are independent of m, and 
the transformation n — t- n — x is admissible. It reduces this case to the case I.b. 

If (yf 7^ 0, there exist two very cumbersome C-integrable equations being the symmetries of 
equations (I2.48P and (I2.49p . respectively. 

In the case II equation (I3.32p possesses at least one generalized conservation law, and thus, 
in accordance with (I3.39p . we can write Ai and A2 as 

59/0 ,5 dfo 



(3.41) 



Ai 



A. 



2/0 dui 2/0 du 

As a result, equation (13.321) casts into the form 



Ml, fo = fo{u,ui). 



(3.42) ut = U5 + - (In /o)xM4 + ^43^3 + (^4^2 + ^5^2 + ^6)^3 + Aju^ + Asul + Aguj + A10U2 + An, 
where Ai = Ai{u, ui). The first canonical conservation law for this equation is trivial, 

2 ax 

Second integrabihty condition (I3.22p is reduced to 



Id 



d 
di 



Pi 



~ ulfo 




2 dA. 



fo fo dui 2/3 \dui 



OA 




Equating the coefficient at m| to zero, we obtain 



(3.43) 



A, 



OA, 



du, 2/2 



dfo_ 
dui 



u\ Zi 



ci/o 



ulZ2 



0(2) ~ 0. 



where Ci is an integration constant. The function Zi is linear in M2, this is why the identity 
Zi = implies two equations, 

2 

+ IOA3/0- 



(3.44) ci 

(3.45) ci 



25/o?4-45 



dui 

25 fo ui 



dfo 



dui 



dui 
dfo 



-30 Ml +5/0— 
OUiOU OUi ou ou 



dui 
(3 + 2^1^3) 



2/0' Ml 



= 0, 

OA, 
du 



3 f^A, 
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The function Z2 is cubic in U2 and therefore the identity Z2 = imphes four equations 
involving also the factor ci. This is why it is natural to consider two cases, Ci = and ci 7^ 0. 
Moreover, in view of Lemma 9, there appears one more fork, Ai — 2A3 = or — 2A3 7^ 0. 
Thus, we have the four cases 



II. a. ci = 0, A3 



Il.b. ci^O, A3 



2 ^' 
2 ' 



II.c. ci = 0, 
Il.d. ci ^ 0, 



A3- 
A3 



2^1 + fu 



where /i = /i 7^ 0. 

Case II. a. In this case the density in condition f l3.23p can be written as 

3 



P2 ~ U2 /o 



5/, 



° duf dui duf 



dfo 
9m 1 



16 



+ P2M. 



By Lemma 9 the coefficient at should vanish, and it determines the function A 



7, 



(3.46) 



A, 



16 /o^ 



d'fo . . dfod'fo fdf, 



"^^ duf ~^ dui du\ 



dui 



In view of f l3.46p . forth integrability condition f l3.24p is reduced to 



d n d 

dtp-' ^ ""'^'Tx 



fo 



duf 



2/0 



dfo 



Pmi + ulu3{P2U2 + P3) + ulO{2) + 0(3) ~ 0, 



where Pi are some functions of first differential order. Equating the term at u\ to zero, we obtain 
/o = {cu\ + a{u)ui + /3(m)) ^ 1 where c is a constant, and a and (3 are arbitrary functions. As 
a result, the equation Pi = holds automatically, and P2 = yields c = 0. Thus, 

5 a , 5 a'u\ + 13'ui 

2 awi + (3 



fo = {a{u)ui + f3{u)) \ Ai 



2 aui + 

In view of transformation formulas fl3.33p for Ai and A2 one can see that the change u — )■ (f{u) 
allows one to simplify /q. If a = 0, without loss of generality we put /o = 1; if /3 = 0, without 
loss of generality we can let a = 1; if aP 7^ 0, we may assume that P = a. 

Thus, there appear the following three non-equivalent cases. 



II.a.l. /o = l; II.a.2. /o 



II.a.3. /o 



a(u) 



Ui Ui + 1 

Case II.a.l. The identities Ci = 0, /o = 1 lead to the relations Ax = A2 = A3 
Aj = 0. Third integrability condition f l3.23p reads as 



A, 



— r-^u'^— (^A - ^ 

dt ^dx\ ^ dui 



+ -ulA^ 



3A« 



dA^ 
dui 



P2{U3) ~ 0. 



In this expression the coefficients at u\ and U3 should be equated to zero. At the same time, 
the density in condition f l3.24p reads as 

dA^ 



Pa 



P2{U2). 



dui duf 

The coefficient at U2 should vanish by Lemma 9. The mentioned three identities imply As 
A8{u), = 3 {As + C2)ui + qi{u); C2 {As + C2) = 0, C2 qi = 0, where C2 is a constant. 
In view of the above results we find 



Pi Ul A'g + P2{U2] 
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By Lemma 9 it yields = 0, and we hence get 

As = C3, A = 3 (C2 + C3)mi + qi{u); C2 (c2 + C3) = 0, C2 qi = 0. 
Now conditions (13.221) and (I3.24p are written as 

r + ^(^2 - 4hi - ^ C3 gi Ml - ^ + Psiu-s) ~ 0. 



c/ 2 d fdAg _ d'^A. 

dt^^ ^dx \ dui du 



Equating the expressions at u\ and u\ to zero, we find Aq and Ag, 

3 3 1 

= {q'l + C4)ul + q2Ui + q3, Ag = -{cj - cl)ul + -C3qiuj + -(C5 + 10 + qf)ui + q^, 

where qi = qi{u) are arbitrary functions. 

Then it follows from the third and fifth integrability conditions that 03 = 02 = 0. The third 
integrability condition determines the function Aiq as a polynomial of third degree in mi, and 
the forth integrability condition determines the function An as a polynomial of fifth degree in 
Ui- In order to determine the coefficients of the polynomials Ag, Ag, Aio, and An we check 
10 integrability conditions. This work, being technically not difficult, requires the examination 
of a great number of options while solving equations. The result is S'-integrable equations 
(13. 2p - fl3.9p . as well as integrable equations being symmetries of the equations fl2.10p - (12.1 3p . 

Case II.a.2. The above obtained formulas for Ai, A2, A3, and also (I3.43P and (I3.46P remain 
to be true. Substituting there ci = and /o = ^, we obtain 

5 , , 5,5, 35 

Ai = -7—, A2 = 0, A3 = -- — , A4 = Aj 



2ui " ' 4ui' ' ' 16^1 

It is easy to check that integrability condition (I3.23P can be written as 
d 



dtp' ^ 



+ 



_i / , OAs dA^ d^AA OAs dA^ d^A^ 

U2U^ [QAs + Qui- h 2ui— +6ni— h 3— 2ni- — — 

\ oui oui ouf J au ou ouiou 

1 o of, dAs od'^As o^Mg (9^5 dA^\ 

+ ulQ{u,ui) + P^iu^) ~ 0. 
Moreover, 

P3 - ulu-,' [a + 12Asu,-2 uf^ - 3 + 4 + P2{u2). 

The coefficients at u\ and u\ U2 in (I3.23P as well as the coefficient at u\ in should vanish. It 
gives us four equations, whose solution reads as 

A , ^2 . gi 2g2 

A5 = gi + As 



where g^ = qi{u). A slightly more cumbersome integrability condition (I3.22p yields 

Aq = C2 + q^ui + 2 c^y/ul + 2 ggn^^^ + q[u\. 
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Due to these results, conservation law fl3.24p becomes 




qj + 15q[ 5g^ + gig2 3 .3 

^-i'^^"' 

2 1 
-y/ul (5 ^2 + ^1^2)' - 7 Wi(2 qiq'i + 15 q'l 




duidu du 5 5 

+ P3M ~o. 

The terms with ul should vanish that implies 

A 1 qI Qa Qa I ^ I — / N C3 3c2 

^9 = C4 - o ^3 + 77: + ^ + — Ml + giMl + 7^VMr(5 92 + ^1^2) 



10 15 ' 25^ ^ ^ ' 4ui 

Then from conditions fl3.22l) - fl3.24p we find A^q and An, but we do not write these expres- 
sions because they are bulky. 

To specify constant coefficients and the structure of the functions qiiu) we check ten integra- 
bility conditions. These conditions are satisfied by the equation fl3.12p and an equation being 
a symmetry of equation (12.181) . 

Case II.a.3. The way of arguing in this case is exactly the same as in II.a.2, but there are 
small differences in the formulas. General in the case II. a formulas become here 

5 5a'ni ^4 ^4 

where a = a{u) is an arbitrary function, ,^ = Ui + 1. 
Integrability condition f l3.23p reads as 



d 2 . ^ f . .dA^ dA^ 15 a' \ ^ f dA^ ^ ,dA 



15 



dt \ dui dui du\ 2aC,'^J \ du du 

o d'A, ^ 2a^dM _ gO^^^ - ^ + (2 a'' - aa")\ + 

duidu a dui a 2 0^^^ / 

+ Qi{u, Ui)u\u2 + Q2{U, Ui)ul + P2(m3) ~ 0. 

This condition together with the formula for the density of conservation law fl3.24p 

and by the relation Qi = lead us to four equations with the solution 

A ,92 . qi 2g2 , 4 a' 

^5 = ?! + ^8 = -777 



where qi = qi{u). Then, by integrability condition f l3.22p we determine the Aq, and by inte- 
grability condition fl3.24p we find the function Ag. After that we determine Aiq and An by 
conditions fl3.22p - fl3.24p . All these expressions involving arbitrary functions of u are rather 
bulky and we omit them. 

It follows from the fifth integrability condition — p4 ~ that a' = 0, qi = 0, q'o = and so 

dt 

forth. Only in the An their remain two arbitrary functions of u. The integrability conditions 
5-7 yield a vast system of algebraic equations for the constants and for the two remaining 
functions. It follows from this system that all the functions Ai are independent of u. This 
is why we can apply a transformation u ^ u — x leading to the case II.a.2. Thus, in the 
considered case there are no new integrable equations. 
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Case II. b differs from the previous ones by that canonical conservation law fl3.22p has the 
second order. It follows from condition f l3.22p that 

5 

/o = —7- — {ul + a{u)ui + ah' = 2a'b, 

and the functions A^, A^, As, and Ag are expressed in terms of /o, 



15 d'fo , 59/0/, df, 



" ju ^ _j_ " u I y "Jo ^ 

2 /o dudui Jq du \ dui 

A ci^/o 5 d^h 35 d^. dh 
^ 4 9mi 8/0 dul 32/2 ^^2 8f§ \du 

A ^ (^Af '.^f^ \ ^'^0 ^ f ^'^0 ^'^o^^o 



24 /q \ 9ui / dudui 12/0 \ duduf dul du 

3 du 24/^ 9m (9mi \ dui 

The formula for is omitted since it is bulky. 

Taking into account formulas fl3.33p and an explicit form of the functions Ai and ^42, one 
can observe that if functions a and h are nonzero, then we can make them constant by a point 
transformation u — )■ ^{u). If a = 0, then up to a point transformation we have either fe = or 
6 = 1. If a 7^ 0, by letting a = 2 we get h' = 0. Thus, there are the following three possible 
cases, 

5 „, _ „ 5 



Il.b.l. f, = -——; II.b.2. / 



2ciM?' ■ ■ ■ 2ci(Mf + l)' 

II.b.3. f, = --^{(u, + iy + c)-\ 
2ci 

where c is a constant. 

Case Il.b.l. It follows from conditions (KT2^ and (K^ that 

^5 = ^8 = 0, ^6 = 02 + qiul + g2M^^ ^7 = ^1^ , 

o 

13 5 

^10 = q'M + g2^^^ ^9 = "2 91^1 - 2 '^2^1^^ ~ 2 ^2^r^ 

.2 \ I "^2 3 3 1 1 2„ -3 



^11 = ^ (^^l + ^1 J + ^ ^1^1 - 5 C2 92^1 ~Y0^^"^ + (^YS ^^^2 - 392J "1 + ?3, 

where gi = C2 is a constant. 

The check of the conditions 6-10 shows that there exist only two integrable equations, 
which are the symmetries of equations (12.151) and (12.171) . 

Case II.b.2. It follows from conditions I K22^ and fl3:Mll that 



2 



2 , 5 19-9^2 



^5 = ^8 = 0, Aq = q + C2Ui\ uj + 1 + {3 q + C3)u^, A-j = -ux , „ , q = q{u) 

V 8 [u\ + 1)'^ 

3 2 + Cs Co 1 \ , , , , / 9 



^9 = 2,1 + / o ~ o'^2\/^i + 1 ~ o (3g + C3)ui, = g'ni(3ni + 2), 

^11 = ^ C2(3 q + C3)(m? + 1)^/' - i C2(2 g + C3)(m? + 1)^/' + ^(3 + 2 C4 g) mi+ 

+ |^(l0g" + (3g + C3)' + c2)nf + ^(5g" + 6g2 + 5c3g + c2)n? + C5, C2g' = 0. 
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The check of the conditions 6-10 shows that there exist only two integrable equations, 
which are the symmetries of equations (12.141) and (12.161) . 

Case II. b. 3. Second integrabihty condition (I3.22p allows us to show that all the functions Ai 
are independent of u. This is why by the transformation u ^ u — x equation (I3.42p is reduced 
to the equations from the cases Il.b.l if c = and II. b. 2 if c 7^ 0. Thus, in the considered 
case there is no new integrable equations. 

Case II. c. The density in (13.231) is equivalent to the cubic in U2 expression (I3.35p . The 
condition tp^ = allows us to express Aj in terms /o and /i. Then we find 

^P2 ~ iZiU2 + Z2) ul + 0(3). 
From Zi = 0, Z2 = we deduce two equations of the form 

= Fiifo, /i), = -^2(/o, /i), 

which can be integrated explicitly. Substituting Aj and Ag in all the expressions, we get 
P3 ~ awl + 0(2). Since 2 A^ — Ai = 2/i 7^ 0, by Lemma 9 we have a = that yields the 
Riccatti equation 

|{^ = ^l(/o)/l' + ^2(/o)/l+^3(/o), 

where (p2 and (fs depend both on /q and on the first and second order derivatives of /o w.r.t. 

Ml. 

As above, we obtain ~ m|((5iM2 + Q2) + 0(2) and consider the equations Qi = 0, Q2 = 0. 
The second of these equations determines A^, and the first implies an ordinary differential 
equation with the derivatives of /o w.r.t. ui containing fi. By the forth integrabihty condition 

~ uliPi U2 + P2) +ulPz + ulu^{Pi U2 + P5) + u\{P^ u\ + Pt U2 + Pg) + ^^3^9 + 

+ m^(Pio ul + Pii ul + P12 U2 + P13) + ^30(2) + 0(2) ~ 

we obtain equations Pj = 0, z = 1, . . . , 13, among those there are many equations involving 
only /i, /o, and the derivatives of /o w.r.t. ui. Expressing all the derivatives of /o from some of 
the equations and substituting them in other equations, we obtain the contradiction /o/i = 0. 

It means that under the conditions II. c there exist no integrable equations. 

Case II. d. We recall that in this case the functions Ai and A2 read as (I3.4ip that ensures the 
triviality of the first canonical conservation law. The function is given by formula (I3.43p . and 
since ci 7^ 0, we have extra two equations (I3.44p and (13.450 . Moreover, A^ = y4i/2 + /i, /i 7^ 0. 

After the exclusion of /I3, equation (13.441) casts into the form 

P ^« If - ^« ( t ) ' - ^ V« t - ^« ^0 1 - «• 

and equation (I3.45P allows us to express A^ in terms of /o and /i, 

= 2;^ 5^ - 3;^ 5^ I' - ' - ' ^0 " 3 

In view of said above, second integrabihty condition (I3.22p becomes 

^pi ~ ul{Zi ul + Z2UI + Z3 U2 + ^4) + Z^ul + Pe{u2) ~ 0. 
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We express Aj from the equation Zi = 0, 
(3-49) A, J-^ S - ^ f 1^ + 84 /o h] 



112 dul 1568 V dui ■""'J duj 
and Ag from the equation Z2 = 0, 

(3.50) As=-/, -i^-^-i^^ + Y^(30c,.,/o-7A)^ 

-168r"^9^-''°^° + ''"^^°^V^"^^'' W ^ 

+ 5 /o ^ /272 — 63)^^^^ — ■^0 ^ ^/l (^31 -^^ ^■^o — 7/c 
504 dui du 63 c^-u \ dui 

By the equations ^3 = and Z4 = we can express the functions and Aio, respectively, 
in terms of the functions /o, /i, and their derivatives. We omit these expressions since they 
are cumbersome. 

The equation Z5 = is an ordinary fifth order differential equation for /o w.r.t. the variable 
Ml. Other implications of the second integrability condition involve the derivatives of the 
functions /o, /i, Ag, and An w.r.t. two variables mq and ui and are too complicated for the 
analysis. 

Then, by Lemma 9 it follows from expressions fl3.35p for p2 and ps that '?/'4 = and ip^ = 
These two equations read as 



(3,51) 70 /„^^ -/„ ^405^^-28 A /.^^ 



(3.52) 25 /o^ - 10^ ^5^ - /o /. J + /^(15 c. /o + 28 ft) ^ 0. 

Expressing the second derivative of /o from fl3.52p and substituting it into (13.511) . in view of 
(I3.47P we obtain /o = — 4/(5ci)/f. After the exclusion of /o, equations (I3.47P and (I3.52p are 
reduced to the equation 

(3.53, 25A^-75(f)%iO/?f ,8,J.o. 

and equation (13.5 ip together with the mentioned equation Z5 = are the implications of 
equation (13.531) . Substituting /i = 5/(4/) into (I3.53p . we obtain the equation 

(3.54) 2A^,^/^ . 9f 



dui \ dui J dui 
whose general integral is written as 

(3.55) (/ + Mi + a)2(2/-Mi-a) + 6 = 0, 

where a and h are arbitrary functions of the variable u. 
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In view of all obtained results including differential consequences of equation f l3.55p . it is easy 
to check that 

^pi ~ 4 r'\3a + 3ni - 5 f){3a'b - ah') + P,{u2), 

where the prime denotes the derivative w.r.t. u. Thus, 3 a' b = ab' that implies a = cb^^^, c = const 
iib^O. 

This result allows one to convert both these functions into constants by a point transformation 
u = f{v). Indeed, since 2 /i = 2 — Ai, in accordance with formulas fl3.33l) . fi{v) = (p'fi{u). 
Therefore, the function / ~ f^^ is transformed by the law f{v) = ip'~ f{u). Making the 
transformation in equation (13.551) . we obtain 

(3.56) [f + vi + a{u)ip'^^]'^[2f -vi-a{u)ip'~^] + b{u)ip'~^ = 0. 
If a = 6 = 0, then no transformation is needed and we have 

(/ + ni)2(2/-ni) = 0. 

If 6 = and a 7^ 0, then letting ip' = a, we reduce equation (I3.55P to 

(/ + ni + 1)2(2/ -ni-l) = 0. 

If b{u) 7^ 0, then a{u) = kb^^^{u), where k is a constant. Choosing if' = b^^^, we obtain 
equation (13.551) in the form 

(3.57) (/ + ui + a)2(2/-ui-a) + l = 0, 
where a is a constant. 

Thus, up to a point transformation, the quantities a and b in (I3.55P are constants, and the 
following three cases are possible, 

Il.d.l. f = -ui- a; II.d.2. / = ^{m + a); II.d.3. f{ui) satisfies fl337p . 

In each of these cases the parameter a takes one of the values, a = or a = 1. 

Employing equation (13.541) . we can exclude higher derivatives of / from the expressions for 
the functions A^ found above. It leads us to rather compact expressions, 

^ = -^/', ^ = 0, As = ^(1-2/0, ^4 = ^(16/'^- 3), 

^5 = 0, Ar = -^f-'i2f'-l){28f" + 20f' + l), As = 0, 
lb 



K2d'A, 1 d'A, IdA, 3f 

2^ ^ + 4^^^^ + - 4 9^ - 2/^^' 



^2 d^Ae , 1^. ,,,,, , dAe 1,^ , , ^,,,,dA, 



Ao = ^ir^rW + 0^1/(2/' + - 7t(/ + ^1 + 2//' 



dudu\ 2 dudui 2 du 

and these formulas hold true for each of three cases Il.d.l, II.d.2 and II.d.3. 

Case Il.d.l. If a = 1, then it follows from the integrability conditions that Ai = Aj(Mi), Vi. 
Thus, the transformation u ^ u — x is admissible, and we arrive at the case a = 0. In the case 
a = numerous forks lead to the only integrable equation (13.111) . 

Case II.d.2. If a = 1, then exactly as in the previous case we arrive at the case a = 0, and 
in the case a = we get equation (13.101) . 
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Case II. d. 3. Let us consider this case in more details. In view of the above results, the 
second and forth integrability conditions read as 

(3.58) ^pi ^ulQi + ulQ2 + ui'Q^ + ulQ^ + 0(1) ~ 0, 

(3.59) ~M4(^l + + ^3(^3 + U2Pi) + 

+ ul{P^ + U2P6 + ulPj + U23P8) + ulP^ + P(,{U2) ~ 0. 

Here the functions Qi and Pj depend on uq and Ui only. For the equivalence of these expressions 
to zero one needs the identities Qi = 0, Pj = for all The conditions Qi = 0, P2 = 0, 
P4 = 0, Ps = 0, and Pg = are homogeneous ordinary differential equations for the functions 
Aq{ui), and u in involved as a parameter. The first two equations are of the fifth order, the 
orders of the others are 6, 7, and 9, respectively. By excluding higher derivatives from the first 
two equations, we arrive at the equation 



+ /(/' + mf - + (1 - 3/') A = 0. 



All remaining equations are its differential consequences. The general solution of the above 
equation is given by 

(3.60) Ae = 7(m)(/ + mi + af + 10ujiu){ui + a)f, 

where 7 and u are arbitrary functions. 

Substituting solution fl3.60p into the equation Q2 = 0, we obtain 

(3.61) 70 aco'iui + a)/^ + a 7' hf + U{ui + a)f + 7{ui + aff + f - ui - a\ =0. 
Calculating the resultant of the polynomials fl3.6ip and (13.571) w.r.t. the variable ui yields 

a [34300 w'2(20 cu' + 97')/^^ + 980 a;'(165 cu'^ + 350 w'^ + 37'^)/^ 



77 (930 wV + 2100 co'^ - y^)f + 7'^(210 cu' + 597')/^ - 7'^ 



0. 



Since u and 7 are functions of u, and / is a non-constant function of Mi, all the coefficients of 
this polynomial should vanish. It implies 

07' = 0, au' = 0. 

Next we consider the equations involving An. These are = 0, Qa, = 0, P5 = 0, and 
Pg = 0. First two of them have the second order, from the first we can express d'^Au/dul, and 
d'^Aii/dudui from the second. By = we can exclude An from Pg = 0, and it implies the 
equation for the function u, 

io" = 6uj\ 

Together with au' = it implies that au = 0. 

Excluding higher derivatives of An from P5 = 0, we obtain the equation 

f)A 

= P,if,unl',co')P,~\f,m,^',u'), 

where Pi and P2 are polynomials in the variables / and ui, such that Pi(/, mi,0,0) = and 
P2(/, Ui, 0, 0) 7^ 0. It means that if a 7^ and 7' = = 0, then Aq and An are independent of 
u. It leads to the fact that all Ai depend only on Ui. Thus, the transformation u ^ u — ax is 
admitted and it eliminates a in equation (13.571) . Therefore, it is sufficient to consider only the 
case a = 0. 
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As a = 0, the equations for An become not very bulky, 




Y=Jf (7" + 4 7^) (8 r + 16 «i f + 8 m'f + 3f + u,) 
+ ly (9 r + 18 f + 9 m'f + 4 / + Ml) 
- ^co^iU f + 23 u, f - 31 MiV' + 9 / - 3 Ml), 



dAu 



du 



40/ 



2P + uiP-ui^f + l 



-2/ 



2 (2 /3 + 3 Ml + Ml V + 1) (2 70;^ + 3 y g;) 

2/3 + Mi/2-Mi2/+l 



Integrating the first equationp we obtain 



^11 =^{8p + 2f-m' + 16 Ml + 8 Ml V') (4 + 7") 

+ 1 7' (18 niV^^ + 36 Ml r + 7 f + 18 f - m^) 

- 4 co\2 f + 3f-mp-8 u,^f + Mi^) + « Ml + /3, 



where a and /3 are arbitrary functions of u. It follows from the second equation that a and (3 
are constants. Using the Galilean transformation, we can assume that a = 0. 

Substituting the expression for dAu/du into the equation = 0, we arrive at extra two 
equations for 7 and u. Finally the system of equations for these functions reads as 



Finding all these functions Ai explicitly and employing equations fl3.62p - fl3.64p . it is easy 
to check completely the integrability conditions 1-4. These conditions lead us to the only 
restriction (]u' = 0, where /3 is the constant involved in An. 

If cu' 0, then /3 = 0. If u' = 0, it follows from (K62^ that a; = 0, and from flHIMD 
that 7 = const. In this case the coefficients of equation (13.421) are independent of m, and the 
transformation u ^ u + pt eliminating the constant /3 in An is admissible. Thus, /3 = for all 
u and 7. 

If a; = 0, letting 7 = 5/i, we obtain equation fl3.13p . 

If a; 7^ 0, it follows from fl3.62p that u' 7^ 0. In this case the order of equation fl3.62p lowers 
and we obtain the equation w'^ = Aoo^ + c coinciding with (13.181) . Since uj' 7^ 0, then (13.641) 
implies 7 + 15 0; 7^ 0, and hence one can express 7' from (I3.64p . It allows us to exclude the 
derivatives of the functions 7 and u from (I3.63p . As a result we obtain the equation 



where c is a constant in (13.181) . 

If 7 = — 30ci;, then it follows from (I3.64p that u = 0, which contradicts the assumption. If 
7 = — Soj, then we get the equation (13.141) . and if 7 = —20a;, then we obtain equation (I3.15p . 

^The method of integration is described in Appendix 2. 



(3.62) 
(3.63) 
(3.64) 



7"' = 87'w + 47a;', 

(7 + 15^)7' + 10(7 + 10uj)uj' = 0, 



(3.65) 



(7 + 30 w) (7 + 5 co) (7 + 20 tu) [(7 + 20 w) (7 + 5 a;)^ + 125 c] = 0, 
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Consider the case 
(3.66) (7 + 20a;)(7 + 5a;)2 + 125c = 0. 

Cubic curve ( I3.66P is rational and is parametrized by the substitution 

CO = u + coj^'^ , 7 = — 5 cw"^ — 20 cD, 

where c = —27 c. Substituting these expressions into fl3.62p - fl3.64p . we find that u satisfies 
equation (13.181) with the constant c instead of c. Substituting the found functions Ai into 
equation (I3.32p . employing the expressions for u and 7, and redenoting cD -> w, c — )■ c, we 
obtain equation (I3.16p . 

3.4. Differential substitutions relating equations in the list. 

As it was noted in Section 2.4, while calculating differential substitutions, it is useful to 
know the orders of the canonical conservation laws. In Table 2 we provide the orders of several 
canonical conservation laws for equations in list (13.21) - (I3.16p . 

Even densities are not indicated in Table 2, since they happened to be trivial p2n ~ 0. For 
equation (I3.12p the orders of pi and pg are provided for the case of generic constants; if p = 0, 
then pi ~ 0, Pq ~ 0. 

The differential substitutions admitted by fifth order S'-integrable equations are shown on 
Figure 2. 

Below we provide the substitutions for the equations with generic constants. 

^(/ + Ml). 



U2 



(IlaSD^dMl): u = ^ 



(EHID^dMD: u = Hf + m 

equation (13.181) satisfied by u equals c 



Inip. At that, uj 
27 



4(^^ 



H — A2 v^, and the constant c in the 



16 



A?A2. 



Table 2. The orders of canonical conservation laws. For zero order conservation laws we 
indicate in the brackets to what the density is equivalent 



Pi 




(13:^ 




(1331) 




(13TP 


(i3:s]) 


(ra 


pi 


0, (~ u) 


0, (~ u) 


0, (~ 0) 


0, (~ 0) 


0, (~ U-) 


1 


1 


1 


P3 


~ 


~ 


1 


~ 


~ 


~ 


~ 


~ 


P5 


1 


1 


2 


2 


2 


3 


3 


3 


P7 


2 


2 


3 


3 


3 


4 


4 


4 


P9 


~ 


~ 


4 


~ 


~ 


~ 


~ 


~ 


pll 


4 


4 


5 


5 


5 


6 


6 


6 


Pi 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 




Pi 


2 


2 


1 


2 


2 


2 


2 




P3 


~ 


~ 


~ 


~ 


~ 


~ 


~ 




P5 


4 


4 


3 


4 


4 


4 


4 




P7 


5 


5 


4 


5 


5 


5 


5 




P9 


~ 


~ 


3 


~ 


~ 


~ 


~ 




Pll 


7 


7 


6 


7 


7 


7 


7 





(I3JI)|) ^ (13:91): u = In Ml, 
(HOlD^dSmi): n = ln(/ 



Pi — A2, P2 — —Xf. 
ui) — ln(2 ^), c = 



(E:IS])^(E3D: S = ln(/ + + - ln(-4u;Ar'), 



A1A2 
_ 1 
~ 4 



A1A2. 
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dSHSD^dMD: ^ = ^- — 

zj u 

dSHD^dMI): u = ln{f + ui) 

1 



U) jUi. 



- In 09. At that, oo 
2 ^ 



-Ai</? + 4A^V9-^ c = -108 AfA^. 



( ISm ^dSTBD: u = --\nu.. /ii = Ai, /i2 
(STI^dMD: M = ^/i. 

(IM])^(EI2D: u = -U2-ul±3 Aie"ni - A^e^" + Aac"". 

( I3.12j) — 7- (l3.6j) : ?i = y/ul— jj,. At that, equation fl3.12p should involve an additional term 5 j/'Ui. 
(H^D^diSD: u = 2u2 



M^iGAse-^^Mi + Aie^" 



A^e- 



U 
U 
U 
U 



■ u 



dMD^dMD 
dMD^dMD 
(SaD^dMD 
(I3S])^(I33]) 

Moreover, there also exist the substitutions for special values of the parameters involved in 
the equations. 

Example 2. (I33])^(IMD: u 



-Ui 
2 Ml 
Ml. 
Ui. 



Ui 



-Aie"±A2e2", AiA 



In each of the cases 



A2 = or Ai = 0, the logarithmic substitution m — Inw or m — > — - Inw leads us to a first 

order linear equation for u. That is, the function u can be expressed in terms of u by one 
quadrature. 

If here we let A2 = 0, Ai — )■ — Af, then we obtain substitution (I3.9h — ;-( j3.6p with A2 = 0. 

Example 3. (13.13^ — 7-( l3.8p : u = ln(/ + ui), A2 = 0, Ai = fi. Here one can also express u 
in terms of u by one quadrature. Indeed, as one can check easily, third degree curve f l3.17p has 
the parametric representation 



Mi = -(2e" + e-2") 



/ = - (e^ 
J 3 V 



-2v\ 



at that, V 



u. 



Thus, we have ^ = \ j (2 e" + e"^") 



dx. 
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Figure 2. The graph of the substitutions for fifth order S'-integrable equations 



Appendix 1. Discrete symmetries of Weierstrass function uj 

Equations fl3.13p - fl3.16p can be written in various ways. We note first that in the paper 
[9] these equations are written in terms of the function R = f -\- ui satisfying the equation 
2R^ - 3Mii?2 + 1=0. 

Moreover, there exist transformations preserving the form of equation (13.181) . Indeed, con- 
sider the functions u and u satisfying the following equations of the form (13.181) . 

(Al.l) Lo'^ = Alo^ + c, 

(A1.2) u'^ = 4u^ + k, 

where ck ^ 0. It is easy to check that the simplest transformations 

\ a — a) , 1 o 

(Ti) u = a k = c=-a; 

a + 2uj 2 

k 

(T2) cu = (!; + —, c=-27k; 



(T3) 00 



C + a/c u' 

2^^ 



map a solution of flA1.2p to that of flAl.ip . 

The transformation Ti is invertible, and a) is expressed in terms of u by the same formula. It 
is also possible to invert the transformation T2, but the problem is reduced to solving a cubic 
equation. To invert the transformation T3 one has to solve the Riccatti equation. We note 
that formula ( IT2P helps to find the parametrization of cubic equation (I3.66p . 

The superpositions of the elementary transformations Tj lead to new transformations pre- 
serving the form of equation (I3.18p . For instance, 

_ _ 3 27a2(I;2 ^ 27 o 

Ti * To : OJ = - a -\ ; — - , k = — , c = a ; 

2 2(2w + a)(c:;-a)2 ' 2' 2 ' 
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k 27 ku^ 
T2*T2: uj = u + — - , c = 729 k; 

u'^ [u^ + kj^ 

( l,\ Z.UJ — I— I — , t \/ Z.I l,~ UJ f"*^ 

T3 * Ti : u 



a{2u} + a)"^ - 3V2a^ u' a" 



^ ^ ca)^ + A/ca)(a)^ — 2 fc) 
^ = 2(^3 +\)^ ■ 

Moreover, Ti * Ti is the identity transformation, and T3 * T3 differs from T3 just by tfie sign 
of the root ^/c. Thus, the, equations fl3.13p - f l3.16p can be written in an infinite number of 
ways different from the first impression. 

Appendix 2. Explicit integration of functions depending on ui and / 

To check the integrabihty conditions of equations (13.131) - fl3.16p we need the table of the 
integrals of rational expression R{ui,f), where the function / is defined by equation fl3.17p . 
These integrals can be found by the rational parametrization 

/A O.N ^ 2w^ + l 

of curve fl3.17p . The parametrization allows one to convert the integral of an irrational function 
R{ui, f) of the variable ui into that of a rational function of the variable w, 

The answer can be written in terms of original variables Ui and / by the formula w = f + Ui 
implied by flMJj) . 

To check the integrabihty conditions we have made use of the integrals 



y"</'"c/Mi, n = 0,1,2; -5 ^ m ^ 11. 



For instance, 

lfdu, = \ul - f, I ulf du, = ^ (8/^ + 14 + ful - u,) , 

jujdm = ^ {2f + fu, + 2 ful - 2 ln(/ + m)) , 

/dui , , „ , f Ui dui „ f Ui dui , , „ , , „ 

-^ = 21n(/ + Mi), j J-j± = 2f + m, j ^^ = 21n(/ + Mi) + 21n/. 

To calculate iterated integrals, one should add to the table also the integrals of logarithms, for 
instance, 

j ln(/ + ui) dui = ui ln(/ + Ui) - y - /. 

In the proof of Theorem 2 we have used around two tens of such formulas. 

To check any of given formulas, it is sufficient to differentiate it, exclude /' = 
lower the degree of U\, if needed, by the identities 

M? = 1 + 3u^f + 2/^ u\ = ui{l + 3uif + 2/3), . . . 

implied by f l3.17p . 
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Appendix 3. On recurrent formulas for canonical densities 

Here we discuss the way of obtaining recurrent formulas like (12.21) and (I3.26p . The original 
idea of this method is contained in the work [SS], where a simple method for deducing recurrent 
formulas for the conservation laws of Lax equations was suggested. In the work [10] this 
approach was applied for the linearization of evolution equations and systems. 

For the sake of completeness of the content, we first describe briefly the essence of Zakharov- 
Shabat method. 

Suppose equation (10. ip has a Lax representation, 

^ = [A, L] <^=^ Ut = Un + F{x, M, Ml, . . . , M„_i), 

where by square brackets we denote the commutator of linear operators. For simplicity we 
assume that A = 74(9^, and L = L{dxi^iiU) are scalar differential operators independent 
of 9f, /X is a spectral parameter, m is a solution to equation (10. ip . 

The Lax representation ensures the compatibility of the linear system 

(A3.1) Li) = 0, i)t = Alp. 

We introduce notations for the logarithmic derivatives of the function ip, 

(In^), = i?, {\n^)t = T. 

It is obvious that the functions R and T are related by the identity 

(A3.2) Rt = T„ 

and Rdx + T dt = dlnip. This is why up to a multiplicative constant we have 
(A3.3) ip = exp(^j Rdx + Tdt 

where the integral in the exponent is a curvilinear integral with variable upper limit {x,t). 
Since ipx = i^R, ipt = "^pT, the operator formulas 

hold true. By these formulas we have 

ip~^L{d.j:, /i, u)ij = L{dx + i?, /i, n), ?/'"^A(9^,, /i, u)i}j = A^d^ + R, /i, u). 
Hence, equation (lA3.ip can be rewritten in terms of the functions R and T, 



(A3.4) L{dx + R,fi,u){l) = 0, 

(A3.5) T = A{dx + R,fx,u){l). 

These two equations are nonlinear in R. Their solutions are often sought as Laurent series in the 
parameter /i. Due to (lA3.2p . the coefficients of these series are the densities of the conservation 
laws. 

Example 4. For Korteweg-de Vries equation ut = u^xx — QuUx the associated linear system 
can be written as 

(A3. 6) ipxx -wp- /i^^ = 0, 

(A3. 7) ipt = 4:ipxxx - Quipx - Suxip. 



Formulas flA3.4p . flA3.5p lead us to the equations for R and T, 
(A3.8) Rx + R^ -u- fi^ = 0, 
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(A3.9) T = 4:{d^ + Rf{R) -6uR-3u^. 

Equation flA3.9p can be simplified by (IA3.8|) that yields 
(A3.10) T = ( V - 2M)i? + M^. 

If we substitute the series 

oo 

(A3.11) R = IJ + J^Pnl^'" 

n=0 

into equation flA3.8p and equate the coefficients at the equal powers of /x to zero, we obtain the 
recurrent formula 



1 f d \ 

(A3. 12) Pn+l = -\u6nQ-^^PiPn~i- -^Pn\ , « = 0,1,2,..., 

where 5no is the Kronecker delta. We note that the scale transformation p„ — )■ p„(— 2)"*^ reduces 
the formula to the form provided in the monograph [19]. Let us write down first elements of 
the sequence pn, 

1 1 1/ 

Po = 0, Pi = P2 = -^Ui, p3 = -(m2-m). 

Next, we substitute series (IA3.1ip into equation flA3.10p to obtain the expansion 

oo 

(A3. 13) T = V + ^ OnP"", 

n=l 

where 

(A3. 14) 9n = 4p„+2 - 2mp„, n > 0. 

Since the parameter p is arbitrary, formula flA3.2p defines an infinite sequence of the conserva- 
tion laws 

(A3.15) 1'^" = ^^"' '^ = 1'2,---, 

To obtain the canonical densities p„, it is sufficient to have (1A3.2I) and one of equations 
flMISD or dMl. 

If we use equation (IA3.8p . we arrive again at recurrent formula (lA3.12p . but we lose formula 
(]A3.14p . The fluxes 6n associated with the densities pn can be found from ( ]A3.15p by inverting 

the total derivative operator — (the algorithm was discussed in Remark 4 on page fTTl) . 

ax 

For further reasoning it is more important to understand how to get the canonical densities 
from equations flA3.9l) and flA3.2p . Since equation flA3.9p does not involve the parameter, we 
introduce the parameter apriori and we can choose the structure of the expansion for R as we 
wish. If, for instance, we assume that R is the Taylor series 



n=0 



where p is the parameter, then 

oo 

T = 5^e„p^ 



n=0 
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where the coefficients On are determined by equation (1A3.9|) . It is easy to check that 

" (f d 

6*™ = 4 ^ Pipjpk - 3 Ui5nQ - 6 UPn + ^'T^Pn + 6^ ^ PiPj, 


where we have used the notations for the sums introduced on the page Since in the left 
and right hand sides of this formula the unknown functions On and pn appear simultaneously, 
it does not help for calculating the conservation laws. 

The situation changes if we postulate the expansion of the function R as the Laurent series 



(A3.16) = /i"^ + ^ p„/i". 

n=0 

In this case by equation (IA3.9P we obtain the expansion for T 

oo 

(A3. 17) T = 4^-=^ + 0^2P~^ + + ^„/i", 

n=0 

and the recurrent formula 

1 . 1 . 

Pn+2 





(A3.18) _i^Jp„,,+ l|;p.p^+liLp„j+_L(6: 



12 3 " ^' - 12 
n + 0_i) 5n-i, 



where n = —2, —1, 0, . . . Let us consider the corresponding series of conservation laws (lA3.15p . 
where n = —2, — 1, 0, 1, 2, . . . . If the conservation laws with the indices i ^ n + 1 are known, 
we find pn+2 by flA3.18l) . and then On+2 by flA3.15l) . and so forth. While finding On+2, we have 

to invert the operator — . Under the assumption that the densities and fiuxes of conservation 

dx 

laws ( ]A3.15[) are explicitly independent of t, this procedure is absolutely algorithmic (see page 
[TTj) . At that, the function On+2 is determined uniquely up to an integration constant. 

The beginning of this recurrence is as follows. According to ( IA3.16p . we have p„2 = p_i = 
0, and this is why by flA3.15p we obtain that the corresponding fiuxes are constant, 0_i = 
12c_i, 0^2 = 12c_2- Then we find po = c~2- Next two densities read as 

1 1 C„2 1 

Pi = 2^ + c-i, P2 = Y^Oo - —u ~ ~ ~ ~ 2c_ic_2- 

To determine ^o? we again have to employ equation ( ]A3.15p as n = that implies Oq = cq. 

It is important to note that the constants Cj appearing in finding the fiuxes Oi are not essential 
since they can be eliminated by the change of the parameter p, 

oo 

(A3.20) p^p + ^kip'. 

i=2 

Consider now an arbitrary evolution equation with one spatial variable 
(A3. 20) Ut = K{x,u,Ux, . . . ,Un), n>l. 

In case (10.11) we have K = Un + F{x, u, Ux, ■ ■ ■ , Un-i)- Denote by the Frechet derivative of 
the function K, 

" dK d' 



dui dx^ 

i=0 * 
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The formal series 

k=—oo 

whose coefficients depend on x, u, u^, ■ ■ ■ that satisfies the equation 
(A3.21) Lt = [K,, L] , 

is called a formal symmetry (formal recurrence operator) of equation flA3.20p . It is known that 
the equation possessing generalized symmetries or conservation laws has a formal symmetry 

Equation (lA3.2ip ensures the compatibility of the following pair of linear equations, 
(A3.23) Lip = Xij, ipt = K^ip, 

where A is a spectral parameter. To this system one can apply the procedure of obtaining 
canonical densities described above. Since the operator L is not known apriori, we employ 
equation flA3.5p . 



Let 



1=0 

oo 

k 



(A3.27) i? = p_i/i-i + J]pfc/i' 

k=0 

then 

n—l oo 

(A3.28) T = /i^" + ^ e.ifi-' + J2 ^k^^''- 

i=l k=0 

Indeed, the minimal degree of fi in the right hand side of identity flA3.25l) is contained in the 
term 

and hence the series for T should begin with the term 6'_n/i~". Since n > 1, then p_„ = and 
^_„=const7^ 0. By scaling the parameter /i we convert into one and obtain (lA3.28p . 

Substituting expansions f lA3.27p . (]A3.28p into f lA3.25p and equating the terms at in the 
equation f lA3.25p . we obtain the first density 



The formulas for several next canonical densities can be found in [9]. 

Let us consider now equation f l0.2p and adduce the deduction of recurrent formula fl2.2p for 
canonical densities following of the above scheme. 

1st step. We write hnearization for equation (10. 2p . 



dV dF f dV dF d dF d 



dx) ~'~ du2 \dx J ~'~ dui dx ~'~ du dt 



ij = o. 



2nd step. By the substitutions 

ip = exp ( I Rdx + T dt \ , where Rt = T^, 
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we obtain the equation with "extended derivatives" 

3 Qrn / \ 2 



d ^ 
dx 



dF 



d ^ 

dx 



dF 



d 



which is equivalent to the relation 

d2 



(A3.29) 



T 



d j^d 

dx"^ dx dx 



dui \dx 



(R) + 



dF 

du-2 



du V dt 



d \ dF ^ OF 

dx J dui du ' 



3rd step. We choose an appropriate expansion for R. The simplest choice is to let 

oo 

(A3.30) i? = /i"^ + ^ p„/i". 

n=0 

Remark 12. Our several attempts to find the expansions with the poles of higher order 

oo 

gave nothing new. If, for instance, we assume for equation fl0.2p R = + Pnf''"', then 



after checking several conditions f lA3.15|) we obtain p2n+i = 0, Vra. It is equivalent to that R is 
expanded w.r.t. the parameter ^ = /i^. Similar results were obtained for some other equations 
and systems as well (see [TT] ) 

Having chosen expansion (lA3.30p . we should accept 



(A3.31) 



n=0 



in order to cancel the terms with in equation flA3.29p . 

For expansion flA3.30p we have p_i = l,p--2 = 0, which implies that and are con- 
stants. Since additive integration constants in the fiuxes are eliminated by the transformation 
of parameter flA3.20p . we let 6'_2 = 6'_i = 0. 

Now, as one can easily make sure, substituting expansions flA3.30l) and flA3.3ip into equation 
(]A3.29I) . we arrive at formula (12.21) with indicated there po and pi. 
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